Gravitation, Experiment and Cosmology* 



ON 
On 

S3 
33 
in 



Thibault Damour 

Institut des Hautes Etudes Scientifiques, 91440 Bures sur Yvette, France 
DARC, CNRS - Observatoire de Paris, 92195 Meudon Cedex, France 



ABSTRACT 

The confrontation between general relativity (and its theoretically most plausible deviations) 
and experimental or observational results is summarized. Some discussion is devoted to the 
various methodologies used in confronting theory and experiment. Both weak- field (solar 
system) and strong-field (binary pulsar) tests are discussed in detail. A special discussion is 
devoted to the cosmology of moduli fields, i.e. scalar fields having only gravitational-strength 
couplings to matter. 
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[-»» 1 Introduction 

During the present century, our knowledge of nature has been drastically deepened by (among other things) the 
discovery of two new interactions (weak and strong) and of a relativistic theory of the gravitational interaction (general 
relativity). At present, much effort is being directed towards uncovering a deeper level of description of nature which 
would hopefully unify all the interactions. A minimum requirement would be to unite the classical description of the 
macroscopic spacetime structure provided by general relativity with the quantum description of the microscopic world. 
From dimensional considerations, it seems guaranteed that any theoretical description encompassing both quantum 
theory (with its characteristic Planck constant K) and Einstein's theory (with its two constants c and G, Newton's 
constant) will look totally different from what we know on the Planck length scale 
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l P = — = 1.616 x 1(T 33 cm , (1) 



corresponding to the energy scale 



E P ={ — \ = 1.221 x 10 19 GeV . (2) 

At present it is difficult to conceive of experiments probing directly physical phenomena at the scale (|]). Except, 
maybe by observing relics of the very early cosmological universe. This line of thought will be discussed below. There 
exists another route which might inform us about the way gravity fits at a more fundamental level within the scheme 
of all interactions. Indeed, nearly all the attempts at unifying gravity with the other interactions predict the existence 
of new long-range, macroscopically coupled interactions appearing as "partners" of gravitation. This is notably the 
case in string theory where gravity always appears accompanied by a scalar field (the dilaton), and possibly by an 
antisymmetric tensor field. [Not to mention the many other fields that arise when compactifying a higher-dimensional 
theory]. Our present theoretical understanding of the generation of mass (i.e. finite range) is much too poor to 
allow one to make any prediction about the range of such possible partners. Therefore, the possibility exists that the 
low-energy effective theory, derived from a more unified theory, contains some other long-range field mediating forces 
between macroscopic bodies. 

In view of this possibility it is important to assess clearly what is experimentally known about gravity, defined 
as the result of all the unscreenable long-range interactions between macroscopic bodies. The present lectures adopt, 
as systematically as possible, a field-theory approach to gravitation and try to summarize what are the present 
experimental constraints on any field-theoretical description of gravity. Beyond giving a catalogue of existing, and 
planned, experiments, we try to extract the maximum theoretical information from present data. The two main 
questions that we address are: 

°* Based on lectures given at Les Houches 1992, SUSY-95 and Corfu 1995; to appear in the proceedings of the 5th Hellenic School of 
Elementary Particle Physics. 



(i) which elements of the present "standard" description of gravity (i.e. general relativity) have been really 
tested, and which have not ? and, 

(ii) what types of new fields with macroscopic couplings could have naturally escaped detection so far, and what 
are the most promising experiments to look for them ? 

To complete the point of view adopted in these lectures the reader is urged to consult the (still relevant) Les 
Houches 1963 lectures of Dicke W, and the specialized book of Will 0. 

In Section 7 below we shall turn our attention to cosmological constraints on scalar fields having gravitational- 
strength couplings to matter. 

2 Methodologies for testing theories 

One can usefully distinguish two complementary approaches for testing the experimental validity of any given theory. 
[Though we will apply the following considerations only to gravitation theories, they have a very general realm of 
validity]. These two approaches can be termed the "phenomenological" one, and the "theory-space" one, respectively. 
More simply, they can be respectively characterized by the two verbs "compare" and "contrast" . Before entering into 
the details of these two approaches it may be useful to view the problem in purely logical terms: Let T denote a 
(scientific) theory, and C some of its (observable) consequences. It is well known that T C is logically equivalent 
to (non C) => (non T). This is the rationale for saying that experiments can "falsify" a theory and the basis of the 
phenomenological approach discussed below. In this approach, experiments have mainly a "negative" value, telling 
us something about a theory only when it is "wrong" . On the other hand, scientists would like to have a rationale 
for saying that they can "verify" a theory. The only logical way of doing so seems to consider the set of all possible 
theories say {T'} and to investigate which subset of {T'}, say {Tc}, implies the same consequences C as T. This is 
the basis of the theory-space approach. This approach gives a more positive value to experiments checking that C 
holds: they tell us that the common features (if any) of {Tc} are "true". 

2.1 Phenomenological approach ("compare") 

Let us assume that we dispose of a general "kinematical" model, containing several free parameters, say {pf hcn °}, 
for describing the structure and evolution of some physical system. By (least-squares) fitting this model to the 
actual observations of the physical system, we can "measure" the values of all the phenomenological parameters: 
pobs = (pf hcno ) best -fit- We can then compare the observed values p° hs to any theoretical prediction concerning the 
parameters pi, as deduced from the current standard theory. The final outcome of this procedure is a set of yes- no 
questions 

pf* = pf^y ? ( 3 ) 

Actually, each observed value p° hs comes out of the fitting process equipped with some error bar, say of bs (corresponding 
to some confidence level, and including both statistical and systematic errors). Therefore, the questions (||) should 
be phrased in probabilistic terms. Moreover, as will be clear from the examples below, the theory never completely 
predicts the numerical values of all the p^s but gives them as functions of some underlying theoretical parameters, 
^theory _ ^theory ^theory y ~\Y nen one disposes of more phenomenological parameters than theory parameters, one can 

eliminate the latter and express the p^s in terms of a subset of them, say p' hcory = /* hcory (p* hoory ), where the index a 
runs only over a subset of the range of the index i. 

An example will clarify the phenomenological approach. In the 18th and 19th centuries several scientists re- 
alized that, independently of Newton's theory, it was always possible to represent the motion of the solar system 
by modelling each planetary motion as a perturbed Keplerian motion, with time-varying Keplerian parameters: a, 
e, lj, i, ft, P, T . Moreover each time-varying parameter in the previous list could be decomposed in secular and 
short-period parts according to p{t) = po + pt + \pt 2 + • • • + S„p W7l cos(w„£ + tp n ). This means that, (nearly) in- 
dependently of any theory, one can represent the motion of the solar system by a list of (constant) parameters, 
say {aJf orcury , a Mcrcmy , w^ ercury , tb Melcur y, . . . , a ^ onus , . . .}. In particular, the fit between that extended Keplerian 
model and the observations yielded a certain value for the "secular periastron advance" of Mercury, say w£^ rcuiy . If 
we now assume a particular theory of gravity, we can (in principle) compute the theoretical value of ^Mercury in terms 
of the other parameters in the model [in the process we must use auxiliary relations to eliminate some not directly 
observable theoretical parameters, like the masses of the planets, in terms of the directly observable parameters of 
the phenomenological approach]. Finally, we can compare ^Morcury 8 thcoiy to cj obs . As was discovered by Le Verrier 
in the middle of the 19th century this comparison exhibits a serious (now more than 90 sigma) disagreement between 
Newton's theory and observations. By contrast the prediction for ^Mercury within general relativity, say w^ rcury , is 
in close agreement (within one sigma) with the observed value. Therefore, one usually concludes that the Mercury- 
perihelion test is invalidating (or falsifying) Newton's theory, but confirming (or verifying) Einstein's theory. The 



problem with this conclusion (besides the fact that it is based on only one test) is that the pure phenomenological 
comparison theory /observations is telling us nothing about which elements of the theory are being tested. Which part 
of the structure of general relativity have we actually checked ? and which parts have played no role in the test and 
have therefore not been probed at all ? Are there other theories which pass also with success the same test ? To 
answer such questions one needs to shift from the phenomenological approach to another one which takes more into 
account the various structures of the considered theories. 

2.2 Theory-space approach ("contrast") 

The idea of this second approach is to embed one's currently preferred theory within a continuous space of alternative 
theories. It is well known that our ability to distinguish color nuances is greatly increased if we bring next to each 
other two different nuances to make a contrast between them. In the same way, past experience has shown that one 
can (sometimes) better unravel the inner structures present in a theory if one contrasts it to a theory which is similar 
but different in some way. To use in practice a space of "contrasting" theories, one needs a way of charting it. In 
the simplest case this will mean that we can continuously label the contrasting theories by means of a finite set of 
real parameters, say {/3 a }- [In more complicated cases the labelling will need an infinite set of real parameters, or a 
parametrization in terms of arbitrary functions]. 

Having, on the one hand, a charted space of contrasting theories (together with the predictions they make) and, 
on the other hand, an actual set of experimental data, we can ask which subset of theories are in better agreement 
with experiment. A standard quantitative criterion for measuring the agreement between a set of data, say {a^ bs } 
together with their estimated one sigma error bars {c° bs }, and a corresponding set of theoretical predictions {a;* l hcory }, 
is to compute the \ 2 ("goodness of fit") statistics. In our case x 2 , for given experimental data, will be a continuous 
function of the (3 a parameters labelling the theories (and therefore their predictions): 



It is useful to imagine the function x 2 {(3a) as defining a hypersurface rising above the finite-dimensional space of 
theories. For instance if there are only two theory labels (/3i,/?2)> the theory-space can be plotted as a horizontal 
two-dimensional plane, say (/3i,/?2) = {x,y), so that z — x 2 {x,y) defines a usual surface in the three-dimensional space 
(x, y, z). The best agreement between observations and theory corresponds to the lowest values of % 2 (which is by 
definition positive). Therefore one is interested in the minima of X 2 (Pa), and their surroundings, i.e. in the hollows of 
the surface z = x 2 ( x > V)- More precisely a convenient way of measuring quantitatively the likelihood for some theories 
to be compatible with the observed data is to consider successive horizontal slices of the x 2 hypersurface above a 
minimum, or equivalently level contours of x 2 m t he space of the parameters (3 a when considering only what happens 
in theory space. To each difference in level above a minimum, say Ax 2 = X 2 ~ Xmin' one can attribute a certain 
confidence level (C.L.), which depends also on the number of fitted parameters, i.e. in our case the number of theory 
parameters that we consider. For example, when there is only one parameter (x 2 (/?i)) the condition A% 2 < 1 defines a 
68 % confidence interval around /3™ m ( "one sigma level" ) and Ax 2 < 4 a 95 % confidence interval ( "two sigma" ) . For 
two parameters ((x 2 (/?i> #2)) the 68 % C.L. corresponds to the two-dimensional region A% 2 < 2.3 in the (3\, fa plane 
and the 95 % one to A% 2 < 6.2. Actually, the convenient link we just described between confidence levels and level 
contours of the specific function x 2 {/3a) obtained by fitting to one particular set of data (the one realized in an actual 
experiment) is a simplification. This simple link holds only in particular cases (e.g. uncorrelated Gaussian noise and 
linear dependence on the /3's), or in the limit of large number of data points. In the general case one should consider 
the best- fit parameters, /3™ ln (those minimizing % 2 ), as random variables inheriting their probabilistic characteristics 
(distribution function in the space of the /3a s) from the ones of the noisy data cc° bs supposed to be an arbitrary sample 
selected from a random process with known statistical characteristics. In other words, a more rigourous analysis of 
the confidence level regions in (3 space would need to use Monte-Carlo methods for generating fake sets of "observed" 
data, and would then study the distribution of the corresponding best-fit f3's. 

Summarizing, the theory-space approach ("contrast") associates to each independent set of experimental data 
some confidence region in theory space at, say, the 90 % C.L.. This immediately raises the following questions: does 
the collection of confidence regions corresponding to the various data sets admit a non-empty intersection ? [If not, 
that would mean either that none of the considered theories is correct, or that some sets of data contain systematic 
errors]. And, if there exists a non-empty intersection what is its shape in theory space, i.e. what are the common 
features of the theories that pass the considered tests ? As we see from the last question, the theory-space approach 
is giving us a handle on what theoretical structures are being actually probed by some sets of observations. 

Let us give an example of the use of the theory-space approach. Eddington introduced in 1923 j| the idea 
that, in the quasi-stationary weak-field context of solar-system experiments, it was possible to chart many possible 




(4) 



relativistic theories of gravitation (different from Einstein's) by means of two (weak-field) theory parameters, (3 and 7. 
[This idea was later extended by Nordtvedt, and Will, (J, [| [| who introduced new weak-field, theory parameters: £, 
ai, CH2, 013, £1, £2, C3j C4- We shall see also below how it has been recently possible to extend the theory-space approach 
to the strong gravitational held regime] . As will be discussed in detail below, the Eddington 7 parameter measures the 
average spin content of the fields mediating the gravitational interaction (i.e. as we shall see the velocity-dependent 
or magneticlike gravitational forces), while (3 parametrizes the cubic vertex of gravitational interaction (3-body force). 
By convention, general relativity corresponds to the values [3 — 7 = 1. 

Let us now reconsider within the theory-space approach the Mercury-perihelion test. A relativistic theory with 
Eddington parameters and 7, say T(/3, 7), predicts the following value for the secular advance of the perihelion of a 
planet with semi-major axis a, eccentricity e, and orbital period P 



■ T(/3,7) _ - Newton 



QttGMq 2 + 2j-[3 



c 2 a(l-e 2 )P 



(5) 



where M@ is the mass of the sun and where (X) Ncwton denotes the Newtonian prediction, which is mainly due to planetary 
perturbations if one separates out the effect of the Earth spin precession, and assumes [to simplify the discussion] that 
the quadrupole moment of the sun is small enough to contribute negligibly. We then see that the comparison between 
the observed, w obs ± <7 obs , and the predicted, to T ^ n \ values of the Mercury perihelion advance defines a certain 
confidence strip in the Eddington theory plane. Present data yield a 68 % C.L. strip approximately given by 

2 + 27-/3 

+ — - = 1.000 ±0.001, (6) 

o 

if one assumes that the adimensionalized quadrupole moment of the sun J2 ~ 2 x 10~ 7 0. By contrast with the 
phenomenological approach which led to a yes-no alternative (in the present case: "yes, general relativity passes the 
test"), the result (|^) of the theory-space approach has a much more precise information content, namely: yes, the 
values f3 = 7 = 1 (obtained in the weak-field limit of general relativity) are compatible with the Mercury-perihelion 
data, but so are all the values of (3 and 7 lying in the infinite strip (||) [e.g. (0 = 5, 7 = 3) or (/3 = — 1, 7 = 0), etc. . . ]. 
Many different relativistic theories of gravitation can pass this test which probes only a particular combination of 
velocity-dependent and nonlinear effects. 

This example exhibits the possibility that the y 2 hypersurface corresponding to a set of experimental data has 
the form of a long, flat valley. This shows the need to perform other experiments to find out where, along this valley, 
stands the correct theory. For example, the experiments concerning the deflection of light by the sun probe the 
parameter 7 independently of (3 and reduce the domain of allowed theories to a small parallelogram around the point 
,9 = 7=1. Actually, the two examples of "classic tests" of general relativity that we just gave are somewhat outdated 
and must be replaced by other tests as we shall discuss below. [The Mercury-perihelion test is inconclusive because 
we have no direct experimental measurement of the quadrupole moment of the sun, and the light deflection test is 
superseded by radio- wave deflection and gravitational time delay tests] . 



3 Testing what ? 

The previous section has exemplified the usefulness of embedding our currently favored standard model of the grav- 
itational interaction, i.e. general relativity, within a continuum of alternative models. The next question that arises 
is: what are the natural extensions of general relativity to consider ? To answer this question we need first to take a 
close look at the structure of general relativity. 

3.1 The two structural elements of general relativity 

Einstein's theory of gravitation rests on two basic postulates: 

i) gravity is mediated only by a long-range symmetric tensor field, g^; 

u) Qjiv couples universally to all other (fermionic and bosonic) fields by replacing everywhere (in kinetic and 
interaction terms) the flat Minkowski metric = diag(— 1, +1, +1, +1) of Special Relativity. [See the Appendix for 
our notation]. 

In technical terms these postulates mean that the total action reads 

Stot = Sgig^} + S m bPm,9nu\ , ( 7 ) 

where the "gravitational" action S g is a functional of g^ v only (without any other long-range field, and without any 
preassigned structure, like f^), and where the "matter" action S m is that of the current standard model of particle 
physics [tpm denoting both the fermionic ( "matter" ) fields and the bosonic ( "interaction" ) ones] in which one replaces 



everywhere the fiat metric (and its associated flat connection) by the curved one g^: — » g^ u , — > V M . [With 
the usual subtlety that one must also introduce a "square root" of g^ v , i.e. a vierbein, for writing down the fermionic 
terms; see the Appendix]. 

The replacement requirement f^ v — > g^ is unambiguous for the (spin 1/2) fermions and the (spin 1) gauge 
fields, but leaves open the possibility of introducing an arbitrary dimensionless parameter in the coupling of scalar 
fields to gravity (^^/gR{g)tf^ Lp). In the case of the Higgs scalar doublet this ambiguity has only unobservably small 
consequences at macroscopic distance scales. 

Let us now turn our attention to the gravitational part of the action, Sglg^]. Weyl J^[ and Cartan H (see also 
Ref. fl(i|| ) proved that the most general form of the action leading to second-order field equations in 4 dimensions was 

SM ^J^gJ^ V~9\R(9) - 2A] . (8) 

The constants appearing in eq. (||) are the velocity of light c, the Newtonian gravitational constant G and the cos- 
mological constant A (with dimension [length] -2 ). Cosmological data indicate that the value of A is at most of a 
cosmological order of magnitude (A < 3(i?o/ c ) 2 where Ho is the present value of the Hubble "constant"). Such a 
small value of A makes its presence unobservable in all non-cosmological gravitational experiments. When discussing 
the latter experiments we shall consider that "general relativity" means eq. (||) with A = 0. 

Another way of justifying eq. (||), with A = 0, as being the unique, consistent description of a long-range 
symmetric tensor field in four dimensions is to follow the approach initiated by Feynman fllf . There is a unique action 
describing the excitations of a massless symmetric tensor field h^ v propagating in a flat, four-dimensional spacetime 
which is irreducible and ghost-free (no negative energy excitations). With a suitable definition of h^ v , this unique 
action reads ^ 

s%[h] = \J^- (v - \ h U») 1° V + d»*h - d afl K - d av h«] , (9) 

where h = h" t O = d^, the indices being raised by the flat metric J M „. Eq. (j|) admits the local gauge invariance 
hf+v — * + <9 M £„ + d^p, the presence of which ensures that only positive-energy excitations propagate. The necessity 
of preserving the existence of a local gauge invariance restricts very much the possibility of coupling to other fields 
and to itself. Work by many authors has shown that there is a unique (modulo field redefinitions) way of coupling 
hjj,u in a consistent fashion |L2|, |l3|, [l4|, [H| [r| [Tt], [In other words there is a unique, consistent deformation of the 
linear gauge invariance of massless spin 2 fields] . This unique answer is equivalent to the expansion in powers of n of 
eq. (0) with g^ v = / M „ + k/i m „ (where k = y / 32TrG/c 4 ). 

It is important to keep in mind in the following the two basic structural elements of general relativity. On the 
one hand, the coupling of gravity to all the fields representing matter and its binding forces in the Standard Model is 
described by a "universal metric coupling" , 

S r n = ^Standard Model 

and, on the other hand, the dynamics of the gravitational field itself (propagation and self-interaction) is described by 
the Einstein-Hilbert action (||). Correspondingly to these two elements of structure, there will be experimental tests 
that probe the way gravity couples to matter (in particular the universal features of that coupling) and tests probing 
the structure and dynamics of the gravitational field itself (e.g. its spin and its range). 

3.2 New, macroscopic fields and couplings 

Alternative theories of gravitation are defined by introducing, besides g^, new, long-range fields mediating the gravita- 
tional interaction. However, one should beware that, apart from tensor-scalar theories, all the "alternative gravitation 
theories" that have been discussed in the specialized literature (notably Ref. Q) suffer from various field-theory 
pathologies: unboundcdncss from below of the energy, negative-energy (ghost) excitations, algebraic inconsistencies 
among the field equations, discontinuities in the degree-of-freedom content, causality problems, etc. The number of 
non-pathological field theories that one can construct is actually very restricted. 

In order to prevent any semantic confusion, let us emphasize that, in these lectures, we mean by "gravity" 
(or "gravitational interaction") the actual, observable interaction between macroscopic bodies which extends over 
macroscopic ranges and cannot be screened by presently known means. With this definition, any field having a 
macroscopic range (say A > 0.1 mm), and coherent couplings to electrically neutral bodies, will be said to participate 
in the gravitational interaction or, for short, to be a gravitational field. The aim of these lectures is to summarize 
what is experimentally known about gravity, and to assess what are the various gravitational fields compatible with 
the present experimental evidence. 



Besides the usual Einsteinian field g^ v , with its nearly uniquely defined coupling to the matter of the Standard 
Model, there is no theoretical shortage of fields that could contribute to mediating gravity. [We consider only Bosonic 
fields; see [jTTf for a discussion of the difficulties arising when using the multiple exchange of massless Fermionic fields 
to generate a gravitational- like force] . 

First, there could be one, or several, scalar fields, say ip , (p 2 , . . . , tp n . Scalar fields can exhibit a rich variety 
of couplings to matter. To quote a few: scalar Yukawa coupling to fermions, gs^pip'tp-i pseudo-scalar Yukawa coupling 
gp(fipj 5 ip, "conformal" or "metric" coupling S m [ip m , A 2 (ip)g^ v ] (which means a universal coupling to the trace of 
the energy-momentum tensor, see below), dilaton-like coupling to gauge fields, <^Tr(_F 2 ), axion-likc coupling to gauge 
fields, <pTr(FF*), etc. Moreover, scalar fields can have an infinite range or a finite one without any restriction on 
their sources (contrarily to massless, gauge fields whose couplings are strongly constrained by algebraic consistency 
requirements, e.g. dvF^ — J M => d^J^ = 0). The existence of at least one sort of scalar (by constrast to pseudo- 
scalar) coupling in the list above suffices to generate a coherent interaction between (unpolarized) macroscopic bodies. 
Many theoretical models have suggested the existence in nature of scalar fields: dimensional reduction, extended 
supersymmctry, dynamical solutions to the strong CP problem, the family problem, or the problem of the appearance 
of particular mass scales, string theory, etc. 

Second, there could exist (one or several) vector fields, B^. If a vector field is massive (i.e. of finite range) its 
couplings to matter are not restricted by any consistency requirement (at least at the classical level). If it is massless, 
or acquires a mass only through spontaneous symmetry breaking, it must be coupled to a conserved current. In 
fact, there are several "unused" conserved currents in the Standard Model and it has been suggested that they could 
correspond to new, macroscopic interactions [|l9| ^0], ^2] . 

Finally, let us mention the possible existence of an antisymmetric tensor field, = —B v{i . This possibility 
was raised long ago by Einstein and others |2^| , and has been revived by string theory |24| . As is expected for any 
gauge field, the couplings of a massless B^ are severely restricted by consistency requirements, thereby disqualifying 
the old "unified" theory of Einstein and its modern avatars [ ^5|. On the other hand, a finite range B^ v can exhibit a 
host of phenomenologically interesting macroscopic couplings |26| . 

Besides the scalar, vector and antisymmetric fields (of any order), there are no bosonic fields which are known 
to have consistent couplings with matter and Einsteinian gravity. For example, there is no known way of coupling 
to gravity, in a fully consistent manner, a second symmetric tensor field, be it massless or massive (see [^, [^8| and 
references therein). 

4 Testing the coupling of matter to an external gravitational field 
4.1 Experimental consequences of universal metric coupling 

We shall refer to the matter-gravity coupling ([Tol), i.e. 



as a "universal metric coupling" : all the different fields ip m entering the Standard Model description of matter and 
its binding forces feel the gravitational influence of the external world only through their coupling to one and the 
same metric tensor g^ u . [We have in mind here the case of a test system, of negligible self-gravity, experiencing 
some background gravitational field described by g^v]. Note that, besides general relativity, there are infinitely many 
examples of gravitation theories exhibiting a universal metric coupling (they define the so-called class of metric theories 
of gravity). Indeed, the tensor g^ v to which matter is coupled needs not satisfy the dynamics derived from the Einstein- 
Hilbert action (||), but could instead be algebraically constructed from other fields having their own propagation 
properties. The simplest example of a non-Einsteinian metric theory is a tensor-scalar theory where the g^ v that 
couples to matter via eq. (O) is of the form 



y/g^R(g*)), and A(tp) an arbitrary function of (p. 

Let us now recall a simple, but useful, general mathematical result about (pseudo-) Riemannian spaces (due 
to Fermi and Cartan). We shall phrase it for the case of a four-dimensional Lorentzian manifold (V, g) of signature 
— h ++: Given any worldlinc C in V (not necessarily a geodesic), there always exist coordinate systems x 11 (i.e. a map 
ip from the abstract V to M 4 ) such that the corresponding components of the metric satisfy all along L the following 
conditions 




(11) 




dg^u(x a ) ' 



. 



diag(-l, +!,+!,+!) , 



(13) 



In other words (using Taylor's formula), the metric components g^(x x ) in those special coordinate systems are equal, 
all over a world tube enclosing C, to the usual (constant) components of a flat metric / M „ modulo terms which are of 
second order in the distance away from £,. It is easy to see that in the particular case where the abstract worldline C is 
assumed to be a geodesic its image in _ZR 4 , L = tp(C), i.e. its coordinate representation in one of the "good" coordinate 
systems x M satisfying eqs. ([l3|), is a straight line (x fJ '(s) = x tl (Q) + su^). In the general case, L is a curved line in J? 4 . 

Let us now consider some physical system evolving according to eq. (^) in some given external gravitational field 
g llVl but isolated from the influence of any other external, non-gravitational field. We assume that the gravitational 
field generated by this system is everywhere negligible. Let us introduce a world tube, say T, of spatial radius ~ d, 
which encloses completely the system at all times. We can apply the theorem ( [l3| ) within the world tube T, i.e. define 
a special coordinate system constructed along some central world line £ (to be identified later as the center-of-mass 
world line of the physical system). In the limit where the physical system becomes very small with respect to the 
characteristic scale of variation of g flu (x x ) we can neglect terms of order d 2 , i.e. we can consider that the spacetime 
metric is flat within T: g^u{x x ) = f^v + 0(d 2 ) ~ /^„. In this approximation, the external gravitational field has 
been effaced within T, and we can conclude that the system will evolve (when viewed in the good coordinates x M ) 
as if it was an isolated system in special relativity. This means in particular that the physical experiments taking 
place within the system will exhibit no preferred directions in space (spatial isotropy) , and no preferred velocity states 
(boost invariance). The local time evolution of the system will depend only on the values of the coupling constants 
and mass scales that enter the usual Standard Model. Moreover, the (special relativistic) center of mass of the system 
will follow a straight line in the x^ coordinates, corresponding to a geodesic in the abstract curved spacetime (V,g). 

Then, by comparing the Lagrangian for time-like geodesies [— m(— g^(x x )x fJ/ x u ) 1 ^ 2 } with the well-known La- 
grangian of a test mass in Newtonian gravity [ ^tox 2 + mU (x, i)] we conclude that, when using space-time coordinates 
[x = ct, x l ) adapted to the Newtonian limit, the time-time component of the metric must be given by 

g 0() {x) = -1 + 2U(x)/c 2 + 0(l/c 4 ) . 

As is well-known, this result allows one to predict that, when intercomparing by means of electromagnetic signals two, 
identically constructed, clocks located at two different positions in a static, external gravitational potential E/(x), one 
should observe a difference in clock rates given by 

l + l[[/( Xl )-[7(x 2 )] . (14) 

In eq. (|lj) the subscript i means "when intercompared" by receiving signals from both clocks at some common location 
which can be the location of either clock or, in fact, an arbitrary point (at rest with respect to the clocks). 

Summarizing, the assumption ( |Tl| ) of universal metric coupling is a very strong one which has the following 
observable consequences for the physics of localized systems embedded in external gravitational fields: 

C\ : Constancy of the constants: the outcome of local non-gravitational experiments depends only on the values of 
the coupling constants and mass scales entering the laws of special relativistic physics. [In particular, the cosmological 
time evolution of the universe at large has no influence on local experiments] . 

Ci\ Local Lorentz invariance: local non-gravitational experiments exhibit no preferred directions in spacetime 
[i.e. neither spacelike ones (isotropy), nor timelike ones (boost invariance)]. 

C3: Universality of free fall: small, non self-gravitating bodies, isolated from non-gravitational external forces, 
follow geodesies of the external spacetime (V,g). [In particular, two test bodies, initially next to each other and at 
rest with respect to each other, fall in the same way in an external gravitational field, independently of their mass and 
composition] . 

C4: Universality of gravitational redshift: when intercompared by means of electromagnetic signals two identically 
constructed clocks exhibit the difference in clock rate (or redshift) given by eq. ([l4|) independently of the nature and 
constitution of the clocks. 

Note that C4 yields a non trivial prediction of universality even when considering two clocks at the same location 
(xi = X2). However, this particular case is contained in the more general consequence C\. 

The consequence C3 is also referred to as the "Weak Equivalence Principle" , and the whole set of consequences 
Cj — C4 is sometimes called the "Einstein Equivalence Principle" j^]. In the present lectures, as we wish to distinguish 
clearly experimental facts from theoretical assumptions we will often shun the use of such expressions. Note that 
the consequences C\ — C4 concern only systems with negligible self gravitational fields. Indeed, the universal metric 
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coupling (11 ) leaves open the possibility for self-gravitating systems to feel the external universe in ways that violate the 
consequences C\ — C4. We give examples of such violations in the following. In fact, such violations are generic among 
metric theories of gravity, and general relativity stands out as one of the very few (probably only two [ p9| ) theories 
for which the consequences C\ — C4 hold true even for self-gravitating systems ( "Strong Equivalence Principle" ) . 



4.2 Non-metric couplings and their observational consequences 

In the previous subsection we discussed the four main observational consequences of the postulate that the matter- 
gravity interaction is described by the universal metric coupling (fll|). This is sufficient for conceiving experiments that 
will test the correct ness of the postulate (^TJ). In other words, we are in position to apply the "phenomenological" 
methodology of Sec. |2.l| above (for instance by comparing the free fall acceleration of test bodies). As was said there, it 
is however useful to go one step further, namely to embed the metric couplings (til) within a larger class of non-metric 



couplings. Indeed, knowing what type of violation of the consequences C\ — C4 of Sec. 4.1 can arise when one changes 
the basic assumption (pd|) helps very much in planning and interpreting experiments. 



4.2.1 Dilaton-like couplings 

Kaluza-Klein theories and string theory naturally introduce couplings between scalar fields and gauge fields of the 
form 
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J d 4 x^ V tv(F 2 ) , (15) 
where the trace is taken over the gauge indices of some Yang-Mills field strength, 

f; v = d^Ai - d u A« + f ahc A\Ai , 

so that tr(F 2 ) denotes g aii g^ v F®gF® v . [The gauge potential A° is geometrically normalized, i.e. it contains the gauge 
coupling constant g\. In a D-dimensional spacetime the metric combination ^/gg a>1 g^ 1 ' appearing in (|l5|) scales as 
\(. D - 4 )/ 2 under a Weyl rescaling g^ v — > Xg^u- Therefore this combination is Weyl-invariant in D = 4. This means that 
one cannot then reabsorb the scalar field <p in eq. (|l5|) by a suitable conformal redefinition of the metric . Therefore 



the term (15) cannot be written in the purely metric form (|l l[) . It is an intrinsically non-metric coupling. 

To investigate the observable consequences of this non- metric coupling term we can note that eq. (p"5|), added 
to the bare kinetic term of the Yang-Mills field — tr(F 2 )/4gQ, (if there is any), means that the Yang-Mills field has a 
field-dependent effective coupling constant, 

-J_=4+/^. (16) 

The result (|l6| ) has two types of consequences: i) it predicts that the locally measured coupling "constants" of the 
Standard Model will depend on space and time (violation of the consequence C\ above) , and ii) it entails that bodies 
of different compositions will fall differently in an external gravitational field (violation of C3). Moreover, the spatial 
dependence of the coupling constants will clearly affect the intercomparison of clocks based on different physical 
principles and will violate the consequence C4 above. To show this let us consider the simple case where the only 
coupling constant which varies is the electromagnetic one a (fine-structure constant). A clock based on a Bohr- like 
atomic transition counts time in units proportional to a~ 2 r e where T e = U/m e c 2 is the "Compton time" associated to 
the electron, while clocks based on fine or hyperfine atomic transitions involve higher powers of oT x . By contrast a 
clock based on the stability of a cavity counts time in units proportional to a multiple of the Bohr radius divided by 
c, i.e. to or x r e , which differs in the exponent of a. 

The amount of spatial and temporal variability of g c ff depends on the mass of the scalar field tp, and on the 
presence or absence of other couplings of ip to matter besides (|l^). For instance, if ip is massless and couples through 



(15) to the SU(3) C Yang-Mills field, it will be generated macroscopically through the fact that the mass of protons 
and neutrons (and thereby that of all nuclei) is currently believed to be mainly made of gluon field energy. The 
cosmological expansion of the universe will then induce a slow time variation of ip, reflected in a corresponding secular 
change of the strong coupling constant. 

The existence of a close link between the constancy of the coupling "constants" and the universality of free fall 
has been pointed out by Dicke (see for further references). To see the necessity of a violation of the universality 
of free fall in presence of a dilaton-like coupling it is sufficient to note that the classical action describing the motion 
of a test particle, say an atom, reads (in units where c = 1) 

S m = - mds , (17) 



where ds = (—g l j, u (x)dx ,J 'dx v ) 1 ' 2 and where m denotes the total mass-energy of the atom. The latter mass-energy 
depends on the effective values of the various gauge coupling constants, say a,i — 3off(i)/^ 7r wnere i = 1, 2, 3 labels the 
gauge groups £7(1), SU(2), SU(3) respectively. The tp-dependence of the oVs entails a corresponding (^-dependence, 
and therefore a spacetime dependence, of m: 



m(x) = m[<p(x)] = m[ai{(p(x)}] 



(18) 



Varying the action (J!?]) yields the equation of motion 
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where L^ CT denote the Christoffel symbols of g^ v , and d v m = dm/dx u = (dm/ dtp)(d(p/ dx v ) the spacetime gradient of 
the mass-energy ([l8|). The right-hand side of eq. (|l9|) gives, for, say, an atom starting from rest in the gravitational 
field of the Earth, the additional term <5a = -Vm/ra = — Vln(m) (beyond the usual g associated with g 00 ) in the 
free fall acceleration of an atom. If we compare the free-fall accelerations of two different atoms, labelled A and B, we 
find the difference 



sla — &B = — V(ln niA — In tub) = — 
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(20) 



Since different atoms have different field contributions to their mass energy we expect the brackets in the right-hand- 
side of eq. (pOh to differ from zero. 



4.2.2 Multi-metric couplings, antisymmetric tensor couplings and local Lorentz invariance 

The previous subsection has exemplified how dilaton-like scalar couplings introduce violations of the consequences 
C\, C3 and C4 discussed in § 4.1. However, scalar couplings introduce (in first approximation) no violations of the 
consequence C2, because the value of a scalar field is Lorentz invariant. One needs to consider non-metric couplings 
involving vectors or tensors to exhibit gravitational violations of local Lorentz invariance |30| , Q] . 

Let us first recall how "isotropy of space" shows up in a simple physical situation. Let us consider Schrodinger's 
equation for an Hydrogen atom, 
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-ip = Etp 



(21) 



In this context, "isotropy of space" means the invariance of eq. (21) under arbitrary rotations around the origin. This 



spherical symmetry comes from the fact that A = S^dij and r = (8ijX l x^) 1 ^ 2 are both expressed in terms of the 
same Euclidean metric 5ij. At a deeper level, the latter property comes from the fact that the kinetic terms of the 
electron field [4>l p d ^ ~ rwi/jip) an d of the electromagnetic field (/ QM f" u 'FapF^) involve the same flat spacetime metric 
(7^7" + 7^7' 1 = 2/ A " / ). This coincidence in the propagation properties of the electron and electromagnetic fields will 
be, by definition, preserved in the case of universal metric coupling (|Tl|). By contrast if, for some reason, the coupling 
to gravity of %j) and F^ introduces two different spacetime metrics (say a "matter" metric g™ v for -0 and a "field" metric 
g^ v for Fpv) then there will be observable violations of the "isotropy of space". In first approximation it is enough to 
consider constant metric coefficients. Let us use coordinates (a la eq. (13)) adapted to the matter metric, i.e. such 
that g™ v = (so that in particular = &y for the spatial components). In these coordinates, the field metric will, 
in general, fail to have the Minkowskian form. In particular, the spatial components of the conformal field metric, say 
g^ LV = —gj^v/goo (which are the only quantities that matter), will be of the general form gfj = 5ij + hij with hij 7^ 0. 

In eq. ( |2l] ) r will be replaced by rp — (gf^x^) 1 ^ 2 . Keeping only the terms linear in hij (assumed to be very 
small) leads to an Hamiltonian of the form Hq + H± where the unperturbed Hamiltonian Hq is that given by eq. (|2l|), 
while the perturbation reads 
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(22) 

The (first-order) shifts in the energy levels of the atom are then obtained by diagonalizing the projection of the operator 
Hi in the subspace spanned by some degenerate eigenstate of Hq. Indeed, the spherical symmetry of Hq implies that 
the unperturbed eigenvalues are exactly degenerate with respect to the magnetic quantum number m (we do not 
consider here the accidental degeneracy of the 1/r potential). The perturbation Hi associated with hij will lift the 
spherical-symmetry degeneracy. This gives an observational handle on the violation of spatial isotropy induced by the 
assumption that the electromagnetic field couples to a different metric than the electron. 

In actual experiments, one considers nuclear energy levels, rather than atomic ones, and experimental situations 
where the spherical symmetry degeneracy has been already lifted, e.g. by interaction with an external magnetic field. 
Generalizing the calculation above leads to energy shifts in the \I,M > state (where I is the nuclear spin, and M its 
projection on the magnetic axis) 
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(23) 



where the indices A, B label the protons in the nucleus and x\ B = x\ — x % B . 

Evidently, only the trace-free part of hij will induce M-dependent shifts. For simplicity, we approximate the 
nuclear-structure matrix elements appearing in the r.h.s. of eq. ( p3| ) in terms of those of the electric quadrupole 
moment of the nucleus Q lJ — Y,Ae.{x 
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The electric quadrupole moment operator can be expressed in terms of the nuclear spin I (and of Q 
maximum eigenvalue) as 



= Qzz, its 
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This yields an explicit expression for the M-dependence of the anisotropic energy shifts 



{Ex) 



I,M 



(Z-l)eQ 
i? 3 



^ h ss 5ij 



B % B 3 



3M 2 



1(1+1) 



1(21-1) 



(25) 



in which B denotes a unit vector in the direction of the external magnetic field (quantization axis). Experimental 
limits on the presence of such terms will be discussed below. 

Let us complete this subsection concerned with possible theoretical origins for terms like eq. (E2fl by mentioning 
how they could be induced by certain couplings between gauge fields and a massive antisymmetric tensor field. Indeed, 
if gravity is mediated in part by a (finite-range) antisymmetric tensor field B M „, it could couple to gauge fields via 
terms of the form 

-|tr[(S^) 2 ] . (26) 

[Note the necessity of considering a massive B^ v \ the gauge invariance of a massless one would forbid an algebraic 
coupling of the form ( p6| ) .] When considering electromagnetism, and separating out the terms quadratic in the electric 
field E l = one finds that eq. (|2(]) is equivalent to having introduced (as we did above in an ad hoc manner) a 
different metric coupled to the electric field: namely gfj = Si 



hij with 
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(27) 



One should note also that the coupling ( p6|) implies not only a violation of the consequence C2 [including evidently 
the local boost invariance, the external B^ introducing preferred spacetime directions] but also of C3: the coupling 
of B to F-field energy will, like the dilaton coupling, introduce a violation of the universality of free fall at some level. 



4.2.3 Other couplings of matter to scalar, vector or tensor fields and their experimental consequences 

In the previous two subsections we selected some specific types of non-metric couplings to exemplify clear cut violations 
of a subset of the consequences C\ — C4. In the present subsection we wish to show by means of examples that most 
couplings one can think of, involving scalar, vector or tensor fields, will entail a violation of at least one of the 
consequences C\ — C4. 

As soon as a scalar field has Yukawa couplings, gs^p^^, to some of the Fermions that constitute ordinary matter 
it will be generated macroscopically by the matter external to the test system we are considering. Then the coupling 
of this external ip to the fermions constituting the test system will violate C3. Indeed, universality of free fall means a 
coupling to the total mass-energy content of test bodies, while the Yukawa interaction we are considering couples to 
a total scalar charge of a composite body of the form 

C s = Y^9s<^> (28) 

i 

where the index i labels the various fermions, and <> the quantum average corresponding to the state of the body. 
It seems clear that no choice of the basic coupling constants g l s will be able to ensure the exact proportionality of Cg 
to the total mass. Indeed, even if one chooses the coupling constants to the quarks and leptons so that the scalar 
charge of individual protons, neutrons and electrons coincides with their respective mass, the presence of nuclear and 
electromagnetic binding energies will prevent Cs to be proportional to the mass for nuclei and atoms. 

Note in passing that a scalar having only pseudoscalar couplings (gpip'ipjsvp, ^pe tiVpa F^Fpa , . . .) would not 
contribute to observable gravity because ordinary matter will not generate macroscopic sources for such a field (one 
would need spin-polarized bodies, time dependent magnetic field configurations,. . . ) 



Let us consider vector fields, B^. Contrary to the case of scalar fields, there is a big difference between massive 
(finite-range) and massless (infinite-range) vector fields. Indeed, massless vector fields admit a gauge invariance which 
restricts very much their possible couplings. They can couple only to a conserved vector source. In other words, they 
are generated by a conserved quantity, such as baryon number, lepton number,. . .Although many more possibilities 
are open in the case of massive vector fields, current theoretical lore favours the case of initially massless (gauge) fields, 
even if they are to acquire a mass via some spontaneous symmetry breaking mechanism. In that case, the source has 
to be a conserved quantity, but the observable range of the vector interaction can be finite. If we consider ordinary, 
electrically neutral matter, it offers only two possible conserved quantities: baryon number B = N + Z and lepton 
number L = Z (here iV denotes the number of neutrons and Z the number of protons or equivalently of electrons). 
[Evidently, more possibilities would be open if we were to consider more exotic types of matter. This possibility should 
be kept in mind when discussing the "gravitational" effects of dark matter]. Then the total vector charge of, say, an 
atom can be written in terms of a coupling constant gv and a mixing angle #5 as 

C v = g v \cos 5 B + sin 8 5 L] . (29) 

Again, there is no way of choosing 6*5 such that Cy becomes proportional to the total mass M of the considered atom. 
The best approximation is obtained by choosing #5 = 0, i.e. coupling to baryon number only. However, in that case 
nuclear binding energy makes for a non proportionality between B and M at the 10 -3 level (when comparing a pair 
of atoms). 

Let us consider an antisymmetric tensor field, say B^ v . Both extended supergravities and string theory naturally 
introduce a massless B^ v as a partner of g^. The gauge invariance (B^ u + d^Xu — d v x^) of such a field restricts very 
much its couplings. However, this leaves the possibility of couplings of the form 

- \fs aX ^J a d x B^ , (30) 

where J a is any macroscopic current, which does not need to be conserved. However, if B^ v stays massless, it has only 
one (scalar) degree of freedom in four dimensions and the interaction (^) amounts to coupling this scalar to d a J a . 
A more interesting case arises when one assumes (without being able to exhibit natural mechanisms for achieving 
it) that the initially massless B^ acquires a non-zero mass. Under this assumption, B^ has the three degrees of 
freedom of a massive vector field, and the interaction ( |30| ) is equivalent to coupling that vector to J a . This offers an 
interesting alternative motivation for the existence of finite-range vector interactions coupled to macroscopic currents 



J a Pq |. Note that the term (26) written in the previous section assumed an initially massive B^ v . This assumption 
opens the possibility of many more interactions with interesting phenomenological consequences of which ( p6| ) is just 
an example. However, there is at present no theoretical motivation for introducing such a fundamentally massive field 
[Not to mention the fact that the non-perturbative consistency of the generalized interactions considered in Ref. p6fl 
has not been proven]. 

Finally, note that a common feature of all the non-metric field interactions considered above (<pF 2 , tpipip, B^tpj^ip, 
B^, v e^ va ° da^^pij)), . . .) is their failure to produce, in the case of a composite body (say an atom), a coupling to 
the total mass-energy M of the body. By contrast, the special relativistic result M — j d 3 xT 00 /c 2 with — 
(2c/ \/g) 5 S ma ,tter /$gfj.v shows why any field B... which enters the matter action only by modifying the spacetime metric 
{guv = g^v[B..] — ► Sg^v = C'"8B...) couples to the total mass of a body. [We consider here bodies initially at rest, 
/ d 3 xT (H — 0, and in stationary equilibrium, so that J d z T % i = by the virial theorem, see below]. Metric coupling is 
the only known way to generate a coupling exactly proportional to M, i.e. one which ensures the universality of free 
fall (Ma = F with F oc M). In other words, the universality of free fall (or weak equivalence principle) plays a leading 
role among the consequences C\ — C4, and deserves to be tested with the utmost precision available. 

4.3 Experimental results on the coupling of matter to an external gravitational field 



The observable consequences of a universal metric coupling listed in Subsection [4.l| above naturally lend themselves 
to high-precision, null tests. 

Many sorts of data (from spectral lines in distant galaxies to measurement of solar-system isotopic abundances) 
have been used to set limits on a possible time variation of the basic coupling constants of the Standard Model (3^] . 



For a recent laboratory test of a possible variation of the fine-structure constant see [33 . The discovery of the "Oklo 
Natural Reactor" , a place in Gabon, Africa where sustained (7 235 fission reactions occurred by themselves two billion 
years ago, gave data that led to tightened limits on many constants |34|, ||. In particular ref. [Q quotes for the time 
variation of the electromagnetic and weak (Fermi) coupling constants 

\a/a\ < 5 x 10- 18 yr _1 . 
\G F /G F \ < 10- 12 yr _1 ■ (31) 



See, however, the global analysis of ref. |35| which leads to more conservative limits: e.g. \a/a\ < 10~ 15 yr -1 . 
Improving on previous (already very precise) results, recent experiments J3(| have obtained extremely tight limits on 
any possible space anisotropy in nuclear energy levels. These experiments look for time-dependent quadrupolar shifts 
of the (Zeeman) energy levels of nuclei with spin > 1/2 (in practice 7 = 3/2). In terms of the expression (^5|) above this 
means essentially putting limits on hij = hij — ^h ss Sij, assuming that hij remains fixed in a locally inertial coordinate 
system, while B % (direction of the magnetic field produced in the laboratory) rotates with the Earth. The best limits 
so obtained are of the impressive order 

1. 



hij ^h ss $ij 



<1CT 27 . (32) 

Note that even if the conformal field metric gj„ = —9uvl9oo introduced above happens to be isotropic in some preferred 
frame (maybe some mean rest frame of the universe), say [<?y] preferred frame = (1 + e)8ij, the hij entering Earth-based 
experiments will have an anisotropic contribution ~ euV /c due to the motion of the Earth (with velocity v l ) with 
respect to the preferred frame. As one expects v/c ~ (both from our Galactic motion and our motion with 

respect to the cosmic microwave background) the excellent limit ( |32| ) yields the still very impressive |e| < 10~ 21 . One 
should however keep in mind the assumption (used in setting the limit ([52])) that the source of anisotropy is external to 
the rotating Earth. It seems to me that (because of the possible presence of ill-calibrated DC effects) the experiments 
performed up to now do not put any interesting limits on an Earth-generated "anisotropy of space" , as e.g. would be 
the case for the term (^6|) if had a finite range < the Earth radius. 

The universality of free fall has been tested by many high-precision experiments (Bessel, Eotvos, Renner, 
Dicke,. . . ). Actually most experiments do not let the test masses fall but compare the forces needed to hold them in 
place when submitted to the gravitational influence of an external source (apparent gravitational forces in an Earth- 
based frame). Most modern experiments have used a torsion balance, i.e. a thin wire holding (in its simplest version) 
a rod at the extremities of which are suspended two different bodies. This apparatus measures the non parallelism 
of the apparent gravitational forces acting on the two bodies. Depending upon the way the experiment is set and/or 
analyzed the results probe various types of violations of the universality of free fall. For instance the Princeton ex- 
periment J3^, [ij looked for effects linked to the apparent motion of the Sun. This means that it was probing only 
fields with range greater than or equal to the distance to the Sun. In 1986, hints of apparent violations of both 
the universality of free fall (in the residuals of Eotvos' experiments) and the inverse-square law (in mine data) were 
presented as evidence for the existence of an intermediate-range (A ^ 100 m) force coupled to baryon number (when 
considering non strange matter) p2f . This suggestion has spurred many new experiments, especially ones testing for 
possible intermediate-range violations of the universality of free fall. See Ref. | 3q| for a review of the experimental 
situation and a detailed assessment of the constraints on the intensity, mixing angle 65 (see eq. (p9|)) and range of 
any new macroscopic force. Let us only quote here a sample of the present experimental constraints on the fractional 
intensity a, with respect to gravity, of a force coupled to baryon number [i.e. a = —gy/4irGu 2 for a vector interaction 



(29) with 65 = and u = 1 atomic mass unit]: 

< 10 3 for A = lm, 

|a|<2xl0~ 6 for A = 1 km , 

|a|<4xl0~ 9 for A > 10000 km . (33) 

If one considers infinite-range interactions, the direct phenomenological limit on a possible differential free-fall 
acceleration between two bodies is at the level [[39| |4Q] 

|a A -a B |/|a|<3x 10- 12 . (34) 

One should also note that the most recent analyses of Lunar Laser Ranging data (see section 5.3 below) find that 
the fractional difference in gravitational acceleration toward the Sun between the (silica-dominated) Moon and the 
(iron-dominated) Earth is <10~ 12 . 

Finally, many experimental tests of the universality of the gravitational rcdshift have been performed. In the 
1960's high-precision experiments, making use of the Mossbauer effect, verified that the gravitational redshift of a 
gamma ray line over a 22m difference in altitude was given by eq. @ with 1 % precision |4l| . Other experiments 
have used spectral lines in the Sun's gravitational field, stable clocks transported on aircraft, rockets, satellites and 
spacecrafts, or have compared Earth-bound clocks with the natural clock defined by the highly stable millisecond 
pulsar PSR 1937+21 (For references see Ref. |^]). Some null redshift experiments [xi = X2 in eq.(|l4|)l have also 
been performed. The most precise test to date of eq. ( |l4| ) achieved a fractional accuracy on the gravitational redshift 
~ 2 x 10~ 4 p2[ . It consisted of flying a hydrogen-maser clock on a rocket to an altitude ~ 10 000 km while continuously 
comparing it to a similar clock on the ground. 



4.4 Theoretical conclusions about the coupling of matter to an external gravitational 
field 

As we have summarized above, the main observable consequences of the postulate (|ll| ) of universal metric coupling have 
been verified with high precision by all existing experiments. Within the presently achieved experimental resolution 
of many dedicated experiments, there are no observational hints of violations of the consequences C\ — C4. On the 
other hand, subsection L2 above has shown, by way of examples, that all non-metric couplings that suggest themselves 
within the present framework of theoretical physics generically lead to violations of one or several of the consequences 
C\ — C4. Therefore the simplest interpretation of the present experimental situation is that the coupling of matter to 
an external gravitational field is exactly of the metric form (11). 

This conclusion should not however be interpreted as being final. Let us indeed examine critically the theoretical 
weight of the tests reviewed in the previous subsection. The most impressive experimental limit is eq. (|3^). However, 
no really natural couplings violating the local isotropy of space have been proposed [we exhibited (^6|) as an example 
of field couplings violating C2, but it is rather ad hoc and assumes a massive antisymmetric tensor field to start 
with]. The second most impressive observational limit is eq. (|34|). However, as written down in the first eq. (|33|), 
the data behind ( |34| ) allow, e.g., for a new field, with range A = lm, coupled to baryon number with strength which 
can be as large as 10 -3 times that of gravity. There exist several models in which factors < 1CP 3 appear naturally. 
For instance the old suggestion [^l| of a vector partner of g M „ coupled to the (PCAC) mass current of the quarks, 
rriuUj^u + mdd^d + • • •, generates a force between macroscopic bodies coupled (approximately) to the combination 
B — 0.17 L with strength [(m u + 2md)/rn,N] 2 ~ 3 x 10 -4 weaker than gravity. Even the existence of infinite-range 
(massless) fields should not be dismissed. The couplings of such fields are tightly constrained by eqs . (|3l|) [where the 
first limit is ~ 5 x 1CP 8 smaller than the Hubble rate] and eq. (|34|). It has been recently pointed out [}13|, [14[ that such 
small coupling strengths might be natural consequences of the cosmological evolution. In particular, ref. |44fl suggests 
that the dilaton (or one of the moduli fields) of string theory might exist in the low-energy world today as a weakly 
coupled massless field entailing very small violations of the consequences C\ — C4. 

In view of these possibilities, it is important to continue improving the precision of the experimental tests of the 
consequences C\ — C4. In particular, let us mention the project of a Satellite Test of the Equivalence Principle Q 
(nicknamed STEP, and considered by ESA, NASA and CNES) which aims at probing the universality of free fall of 
pairs of test masses orbiting the Earth at the impressive level 5a/ a ~ 10 -17 . Let us also note that there are plans for 
flying very stable clocks near the Sun; the aim being to improve the testing of the gravitational redshift down to the 



10 6 fractional level, i.e. the level where second-order effects cx (U sun /c 2 ) 2 enter eq. ( |l4| ) (Vessot). See Ref. [46 for a 
recent survey of these, and other, projects in experimental gravity. 



5 Testing the newtonian and post-Newtonian limits of metric theories of gravity 
5.1 What are the most natural metric alternatives to Einstein's theory ? 

In the rest of these lectures we shall adopt the provisional conclusion of the previous section, namely that gravity couples 
to matter in the purely metric way jll]). This conclusion seems to leave open many possibilities for alternative, non 
Einsteinian, theories of gravity. Indeed, the physical metric tensor g pl , through which matter interacts with external 
gravity can still be an arbitrary function of many different fields 

g^v = g^lu, V, B^, B^u, ...} , (35) 

for instance 

9llv = A 2 (v)[g% + ai B^B v + a 2 g;„gP a B p B a + a 3 g^B w B^ + ■■■}. (36) 

[We do not include field derivatives in eq. (|3^), e.g. d^ipdvip, because they induce serious causality problems.] However, 
the appearance, besides a basic tensor field g* and one or several scalar fields ip, of vector fields antisymmetric tensor 
fields, etc. . . in eq. ( |35| ) seems theoretically improbable for the following reasons. First, the non gauge-invariance of 
the combinations B^B V and B pa B vp implies (if one wishes to avoid the presence of negative-energy excitations) that 
the fields B p and B pv must have from the beginning a non-zero mass (or finite range), i.e. more precisely that their 
kinetic terms must be of the form — \F 2 V — ^m 2 B 2 and —j^H 2 ^ — j^B 2 ^ respectively, where F^ v = d p B u — d u B pi 
H\fiu = d\B^ v +d ll B u \ + d u Bx f _ t - [It is well known that all other forms for the kinetic terms lead to ghost excitations]. 
This would mean the presence at a fundamental level of the theory of a particular length scale A = 1/m. It does not seem 
very plausible that such a fundamental length scale happens to be of macroscopic magnitude, as is necessary for it to be 
relevant to the topic of these lectures [if A ~ £p, eq. (|lj), there will be no observable consequences of the presence of such 
fields] . A second reason which does not favor the existence of fields B^ and B pv coupling to matter only through eq. (|35| ) 
is that such fields would exhibit no linear couplings to matter. Their source <5<S ma tter/<5i? = {5S m / 5g pv ){8g lljl , / 8B) is 
(at least) linear in B, so that an everywhere vanishing B field is an exact solution of the B field equations. The only 



way the B fields can couple to local matter is through the presence of a cosmological B background, generated by 
putting suitable boundary conditions at the Big Bang. Thirdly, the quadratic couplings of the B fields to matter, e.g. 
a\T 111 ' B^By modify the mass terms in the action, and it remains to be proven that these modifications preserve the 
consistency of the theory. Finally, though the ellipsis in eq. (35) could stand for other types of tensors (like a second 
symmetric tensor field), we have seen above that it seems very difficult to introduce such fields in a consistent way [i.e. 
free of algebraic inconsistencies, discontinuities in the degree-of-freedom content, causality problems, negative-energy 
excitations, etc. . .]. 

In conclusion, the most natural metric theories of gravity are expected to contain only one symmetric tensor 
field, <7*„, and one or several (massive or massless) scalar fields, (p a , a — 1, . . . n, and to couple to the Standard Model 
of matter via a physical metric of the form 

g^=A i { V 1 ,..., V n )g^ (37) 

where A((p a ) is some arbitrary coupling function. Note that if we require from the beginning to have only massless 
fields the drastic consistency constraints on the couplings of gauge fields (see e.g. Refs. [O, [0| [k| |l4|, |l5|, [l6|, [l?], |l^| 
and j2^, |2(| ^7], |2^]) force one to consider only tensor- multi-scalar theories. [We do not consider here massless Fermion 
fields]. 

Finally, let us give a simple physical argument (which is not really independent of the consistency ones given 
above) which shows clearly why tensor-scalar theories are preferred when one assumes that consequences C\ — C4. are 
exactly satisfied. In fact, let us start only from C3, the universality of free fall. The usual reasoning (F^ = Ma&a 
with = g independent from which body A is considered) shows that gravity couples exactly to mass (Fa oc Ma)- 
On the other hand, Special Relativity tells us that the mass of a body (in stationary inner equilibrium) can be written 
either as 

M= i J d 3 xT 00 , (38) 

or 

M = ^ J d 3 x (T 00 - T ss ) , (39) 

where T^ v denotes the total stress-energy tensor (including matter and field contributions). Indeed, the second form 
is a consequence of the virial theorem 

d 3 xT lj = -dl J d 3 xT 00 x l x 1 (40) 

(which follows directly from the conservation laws dJI^ — 0). Eq. ( [icf ) shows that the integrated stresses, J d 3 xT l i , 
and in particular their trace, vanish for a body in stationary state. 

Now, the first form ( |38| ) of the "gravitational charge" suggests a coupling to a massless spin-2 field, h^T^ (the 
consistency of the coupling being ensured by the fact that T^ v is conserved), while the form (39) suggests a scalar 
coupling ipT^. At the linearized level, we thereby expect an interaction of the general form 

interaction = \ J ^ (V + 2a a ip a f^ u )T^ , (41) 

where denotes as above the flat metric and where the index a labelling the various possible scalar fields is summed 
over. Remembering that T^ v is the functional derivative of the matter action with respect to the metric, 

5S m = ^Jd A x^T^8 gtlv , (42) 

we conclude that, at the linearized level, eq. (^l|) is telling us that the gravitational couplings of matter is described 
by replacing in the matter action the flat metric by 

^linearized = f/u) + ^ + 2aa ^a f ^ „ ^ + ^ftf^ + h ^ , (43 ) 

where h^ v is a massless spin-2 field and {<p a } a collection of (massless or massive) scalar fields. The result ( |43| ) is 
nothing but the linearized version of eq. ([57]). The coupling coefficients a a measuring the relative weight of scalars with 
respect to the spin-2 field in the linearized gravitational interaction are just the logarithmic gradients of the coupling 
function A(ip a ) of eq. @, 

After differentiation, the r.h.s. of eq. (Q) is to be evaluated at the background (or vacuum-expectation) values of 
the scalar fields (e.g. <pg = 0), while the constant conformal factor A(ipo) must be transformed away by rescaling the 
coordinates. [More about this below]. 



5.2 The Newtonian limit of tensor-multi-scalar theories and its experimental tests 

Let us commence by denning in full detail the most general class of tensor-multi-scalar theories. The total action reads 



St, 



s n 



+ s v + s 



(45) 



with 



Sn 



d 4 x 



V9* 



4ttG, 



f d 4 x 
J 



■B{<p°) 



(46) 
(47) 



and 



S m — ^Standard Model [V'm, 9 pis] , (48) 

in which the physical (or "Jordan-Fierz") metric g^ u directly coupled to matter is related to the "Einstein" one g*^ 
appearing in the Einstein-Hilbert action ( f46|) (where i?* denotes the Ricci scalar of g* ) by a scalar-field dependent 
conformal factor 

9tut = A\^)g% . (49) 

The universal coupling of matter to g^y means that (non-gravitational) laboratory rods and clocks measure this metric. 
[It would take a purely gravitational clock, e.g. that defined by the orbital motion of two black holes, to measure 
the metric <?*„]. The action ([l||) contains one dimensionful constant G» ("bare" Newtonian constant) and several free 
functions: the n(n — l)/2 arbitrary functions J a b(ip) entering a general (cr-model) metric in the n-dimensional space of 
scalar fields (da 2 = r y ao (p c )dp a dp b ) and the two functions A(ip a ) and B(<p a ) which give the coupling of the scalars to 
the matter, and the self-couplings (potential) of the scalars respectively. The original theory of Jordan-Fierz-Brans- 
Dicke ^7) has only one scalar field and one free parameter ,a. This theory is defined by the choices A[tp) — exp(aip), 
B(ip) = 0, da 2 = (d<p) 2 . The coupling parameter a — din A/ dip (which is a constant in this theory) is related to the 
often quoted parameter lu through a 2 = (2u> + 3) _1 . 

The gravitational field equations corresponding to the action (E5h read 



r; v = 2 lab (<p)d tl <p a d v <p b + 2B(p) g ;„ 



2q* [ T t 
dB 



, 9uv 



In eqs. (BOh, (p|) we have used the notation 



q* = 

rpflU — 

* 

a a (p) = 



4irG*/c 4 , 
2c SS m [ip m ,A 2 g* 



fan 



19* 



dlnA{p) ^ A _ x dA 



dp a 



dip a 



(50) 
(51) 

(52) 
(53) 

(54) 



Moreover, D gt = gt V* V* denotes the g*-covariant Laplacian, 7 the inverse of 7 a &, 7^, the Christoffel coefficients 



of jab, and the various indices are moved by their corresponding metric: T*„ = g*g*gT"^, a a = 7 afc ab, etc 
Note that the "Einstein-conformal-frame" stress-energy tensor ( |53| ) is related through 

T * u = A 6 T^ , s/g^T"^ = JgT" 

(in which T" = g^ a T av ) to the physical ( " Jordan-Fierz-frame" ) stress-energy tensor 

= — SS mbPm,9f^u] 



The latter tensor satisfies 



y/9 $9n» 



W V T^ V = , 



(55) 



(56) 



(57) 



with respect to the g-covariant derivative V M , while it is only the sum of Tt" and of the stress-energy tensor of the 
scalar fields which is <?*-covariantly conserved. 



One sees from eq. ( pQ) that the scalar potential B(tp) introduces an effective cosmological constant in the tensorial 
field equations, A = 2 < B(ip) > . There are very tight constraints on the value of A in ordinary units. One should 
therefore restrict oneself to considering potentials B(tp) that tend to dynamically drive the scalar fields toward values 
(fig at which B(ipo) = 0. From eq. j5l] ) we see that this means, if jab is positive definite, that B(tp) should have a zero 
minimum value at (p a = tpft. Setting by convention = 0, we can then easily write down the linearized approximation 
of the field equations: 



2\R 



* n linearized 



(□* - m* 2 )y a 



(J, V 

q*a a (0)T, 
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(58) 
(59) 



h* 



Here we have expanded the Einstein metric as g*^ u — f* 
Minkowskian form when using some Einstein- frame coordinates a;*) 
field coordinates tp a that diagonalize the scalar mass matrix, i.e 



where /* is a flat metric (which takes the usual 



d* denotes d/dx'i, □* = f*"dt v , and we used 



,*2r6 



Inserting the solutions of eqs. (58) and (|5 
the gravitational fields h^ v and tp a , namely 



the second-order gradients of B(tp) around zero, 
into the action (|4l]) giving the interaction between the matter and 



(h*Tr + 2a a (0)<p a T* 



jives [using the harmonic gauge d"(h* - \b*£ ' f^ v ) = 0] 



(60) 
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L J *l0C U * \ L1 tlV 1 * JlM/J 
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a=l 



a a (0)a°(0)Ti° c (D* - mf )- 1 T, ext 



(61) 



where T* l ° c is the energy distribution of a local system which is gravitationally interacting with the external energy dis- 



tribution T*^ xt . Eq. (61) shows clearly that the metric g*^ v mediates a usual, Einstein-type massless spin-2 interaction, 



while each scalar field mediates a, possibly massive, spin-0 interaction. 

The Newtonian limit of eq. (^lj) consists in neglecting all velocity dependent terms, which amounts to neglecting 
the components T° l and Tl 3 with respect to the time-time components T°°. [Indeed, for ordinary materials T°° ~ pc 2 , 
\T^\ ~ pev, \Tl° | ~ pv 2 where v is some (internal or orbital) velocity]. This yields the following interaction Lagrangian 
between two (point-like) bodies 



L 



Ncwto 



G*A(0) 



2 MiM 2 
ri2 



a=l 



* o (0)a°(0), 



(62) 



where the factor A(0) 2 comes from having rescaled both the stress-energy tensor and the coordinates when passing 
from the Einstein frame x*i to the physical frame x^, such that ds 2 — g^dx^dx" — A 2 ((p)g*, 1 dx*dx" tends to the 
usual Minkowski metric / M „c 



,dx^dx u at infinity. The scaling transformations are x M 



v4(0)x* and eq. 



<EB) must be 



modified by taking into account the coordinate change x*i — ► x^. [As it stands eq. (|55| ) assumes the use of the same 
coordinates in the Einstein and Jordan-Fierz frames]. The quantities appearing in the final eq. ([32|) are all expressed 
in physical units [e.g. M = J d 3 xT 00 /c 2 using physical coordinates x^ and the stress tensor (pq)]. 

The most evident experimental consequence of the result (|62| ) [besides its pure dependence on the total mass- 
energies (equivalence principle)] is the possible presence of Yukawa-type modifications of the usual 1 jr potential. Many 
experiments have set tight constraints on such modifications. Here is a sample of some recent results |p8, 46 , assuming 
the presence of only one Yukawa term (with range A = 1 /m) 



a(0)\ 2 
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10" 4 if 


A ~ 1 cm , 




a(0)\ 2 
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5 x 10~ 4 


if A ~ 1 m , 




a(0)\ 2 
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10~ 3 if 


10 m < A < 10 km , 




a(0)| 2 
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10~ 8 if 


10 4 km < A < 10 5 km . 


(63) 



[Beware that the coefficient of the Yukawa term, here denoted [a(0)] 2 , because it appeared as the square of the coupling 
constant a(0) of the scalar field, is usually denoted a]. For the same reasons that we evoked above in the case of 
composition-dependent interactions, it seems desirable to continue performing experiments, both in the lOm-lOkm 
window, where the limits are not very stringent and in the A < 1 mm window which is very poorly constrained ffSf . 

In the rest of these lectures, we shall concentrate on the case where there are only long-range scalar fields 
(m a — 0). In that case eq. ( |62| ) predicts a 1/r potential between two masses with an effective Newtonian constant 
given by 

G = Gdfl) = G.[A(<po)] 2 [l + a 2 (^ )] , (64) 

where 

a 2 ^) ee 7 ^K(^M ee r b ^ ^ (65) 

denotes the fractional contribution of all the scalar fields to the 1/r interaction. We have made explicit in eq. (^34|) the 
dependence of the effective Newtonian constant (as it can be measured by a local Cavendish experiment) upon the 
background value (or VEV) of the scalar fields. Indeed, as we are now considering the case of massless scalars there is 
no need to assume any non-zero potential function B(ip) for the scalars. In that case, there is no longer a force term 
cx —dB/dip a in eq. ( j5Tj ) driving the scalars to a particular VEV On the contrary, the long-range coupling of the 
scalar fields to the universe at large, now exhibited by eq. (^l]), makes us expect that localized gravitational systems 



will be embedded in a cosmologically evolving background: (/?oW- Therefore we expect from eq. (64) that the locally 
measured gravitational constant will evolve on a Hubble time scale 

(dG/dt)/G ~ Hq . (66) 

Various types of observational data (including binary pulsar data f49|| ) can be used to look for a possible time- variation 
of the Newtonian coupling constant. Let us only quote here the result obtained by a recent re-analysis of the Viking 
data (radar ranging between the Earth and Mars) 

G/G = (-0.2 ± 1.0) x lO^V -1 ■ (67) 



As H — /175 x 75 km/sMpc = h?s x 7.67 x 10~ n yr -1 with — 1 ± 0.33, we see by comparing (66) and ( |67j ) 
that the present observational results are not putting a very strong constraint on the possible existence of a long-range 
coupling to the universe at large. In fact the analysis of the post-Newtonian effects in the dynamics of the solar system 
(see below) are putting much more severe constraints on the existence of extra long-range fields than the present G 
observations. [This is the case if one assumes a universal metric coupling. If, on the other hand, the matter driving 
the cosmological expansion is a new type of (dark) matter which couples differently to a postulated long-range scalar 
field the G observations may provide a significant constraint on the scalar coupling of this dark matter |5C|].] 

To conclude this section devoted to the Newtonian limit let us recall the shameful fact that Newton's gravitational 
constant is one of the least precisely measured fundamental constant of physics. By contrast to Ti, a = e 2 /Tic, the 
particle masses,. . .which are known with a part in a million precision (or better), G is only known with a precision 
- 1.3 x 10~ 4 

G obs = [6.67259 ± 0.00085] x l(T 8 cm 3 g- 1 s- 2 . (68) 

This lack of precision could become very annoying if ever theoretical physics allows us, one day, to predict the value 
of G in terms of other physical constants. Landau [j52| entertained this hope long ago and conjectured that the very 
small dimensionless quantity Gm 2 /Tic ~ 10 -40 , where m is a typical particle mass, might be connected with the fine- 
structure constant a = [137.0359895(61)] _1 by a formula of the type Aexp(— B/a), where A and B are numbers of 
order unity. Recently, 't Hooft f33|| resurrected this idea in the context of instanton physics, where such exponentially 
small factors appear naturally. He went further in suggesting (for fun) specific values for A and B in the case where 
771 is the electron mass. Actually, the final formula he proposed is in significant disagreement with the observed value 



(68). However, keeping his (instanton-motivated) value for B, namely B = ir/4, but taking for A the value (7w) 2 /5 
one can (still for fun) define a simple-looking "theoretical" value for G by G theory m 2 /Tic = (7ir) 2 /5 exp(— ir/4a). 
Using the central values of the 1986 adjustment of the fundamental physical constants fjlfl , this formula "predicts" 
G thoory = [6.6723458 • • •] x lO^cn^g-V 2 , which is in good agreement with the observed value @: G ob 7G thcory = 
1.00004 ± 0.00013 ! Let this exercise serve as a reminder of the potential importance of improving the precision of the 
measurement of G. 



5.3 The post-Newtonian limit of tensor-multi-scalar theories and its experimental 
tests. 

The term "post-Newtonian" refers to the terms in the Lagrangian describing the motion of gravitationally interacting 
bodies which contain a factor 1/c 2 with respect to the "Newtonian" terms ( p2[ ) . There are two types of post-Newtonian 
terms: those which are smaller than ( |62| ) by a factor (v/c) 2 ["velocity-dependent terms"], and those which are smaller 
by a factor GM/rc 2 ["non-linear terms"]. 

The velocity-dependent terms (also called "gravitomagnetic" terms) can be directly deduced from the linearized- 
order action (|6l"| ) by inserting the point-mass approximation of T^ u = J mu^u" 5(x — z(s))ds (after the needed scaling 
transformations). The latter equation shows very clearly that the exchange of massless scalar fields introduces a 
different velocity dependence oc jj dsids2miiri2G(zi — z 2 ) than the one due to the exchange of a massless spin-2 field 
cx JJ dsids 2 mim 2 ui M Mi I/ (2M2M2+/ Miy )G(z 1 — z 2 ) [here G(x) denotes the Green function of the flat-space d'Alembertian, 
and Ui, u% the four- velocities of the two considered mass points]. Similarly to what happened in eqs. (^32|) and ( |64| ) 
one sees from eq. (|6l|) that (in the massless case) the factor a 2 of eq. ( |65| ) will weigh the contribution of the scalars to 
the velocity-dependent terms [see Sec. 3 of Ref. Q for details]. 

Let us now turn our attention to the non-linear post-Newtonian terms oc GM/c 2 r. There are two types of such 
terms. The first type can be easily understood from our previous results. Indeed, eq. (^) showed that the value of 
the gravitational coupling constant measured in a local Cavendish experiment depends upon the ambient values of the 
externally generated scalar fields at the location where the experiment is performed. [The calculation behind eq. ( |64| ) 
considered a gravitating system put in a constant scalar background ipg (of cosmological origin) . Because of the long 
range of the scalars, the scalar background experienced by one body member of an iV-body system is obtained by adding 
the effects of the N — 1 other bodies onto the cosmological background] . Therefore the effective gravitational constant 
ruling the self-gravity of a particular body (say a planet) will be space dependent: G(x) = G(ip(xj), where G(ip) is 
given by eq. (p4|). Now, the total mass-energy of a self-gravitating body depends upon G because of the gravitational 
binding energy, say E grav = Gd( m tot c 2 )/dG ^ 0. The space-dependence of G(x) induces a space-dependence of the 
mass to. As was discussed in Sec. 4.2.1 above, [eqs. (19) and (^p|)], this causes a supplementary term in the acceleration 
of the body, namely 
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(69) 



Such a term is absent in pure general relativity where the gravitational influence of the external universe can be locally 
effaced by introducing Fermi-Cartan coordinates, eq. (|ll|). [See Ref. Q for a general discussion of the "effacement" 
properties present in general relativity, and for references]. The presence in tensor-scalar theories of an anomalous 
contribution to the gravitational acceleration of a body proportional to _E grav /TO,c 2 was discovered by Nordtvedt Q. 
[The possibility of such an effect was first noticed, via the reasoning behind eq. (|69|), by Dicke J55[ ]. 

The second type of non-linear terms are the genuine 3-body interaction terms in the action for gravitating bodies. 
To obtain them one needs to go beyond the linearized theory written down in eqs. (58), ( [59| ) above, and study the 
quadratically non- linear terms in the field equations, i.e. the cubic terms in the field action (45) [Fortunately, it 
suffices to study these terms in the slow-motion limit] . An elegant way of dealing with these quadratic nonlincaritics 
has been recently found both in general relativity |5(| and in tensor-multi-scalar theories [^9|. Let us quote the final 
result for the Lagrangian describing, within the first post-Newtonian approximation, the gravitational dynamics of N 
(self- gravitating) bodies [with positions z A (t) and velocities v A (t); A = 1 — N] 
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where 
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(73) 



with tab = \ Z A ~ z b\ and n AB — ( Z A ~ z b)/ r AB- The physical metric corresponding to the post-Newtonian level of 
accuracy can be written as [using the short-hand notation 0(n) = 0(c~ n )] 
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in terms of the following scalar and vector potentials (□ = f^d^) 
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(77) 
(78) 



[One should keep in mind that the post-Newtonian limit is a combined weak-field, slow-motion expansion, so that 
the error terms 0(n) = 0(c~ n ) in cqs. (f74|) - (76) contain both velocity-dependent terms (or time derivatives), and 
higher-order nonlinear terms]. Besides the (dimensionful) constant G, eq. (|64|), there enters only two (dimensionless) 
parameters in the post-Newtonian limit of tensor-multi-scalar theories: 7 and (3. [They coincide with the parameters 
introduced by Eddington long ago when considering the simpler model of test particles moving in the field of one 
central, massive body ||]. The post-Newtonian limit of general relativity is obtained when 7 = 1 = (3. [Note the 
simplifications of the non linear structure that arise in this limit where 4/? — 7 — 3 = = 3 7 — 2(3 — 1] . 

The quantity 7 — 1 parametrizes the possible presence of non-general-relativistic velocity-dependent terms (see 



eq. (|72|)). From our discussion above it is clear that 7 — 1 must be proportional to a , eq. (p5|), which measures the 
admixture of the scalars in the two-body interaction. More precisely, one finds 



7-1 



1 + a 2 



(79) 



The result ( fZQj ) can be formally generalized to the case where the gravitational interaction is mediated not only by 
(massless) spin-2 and spin-0 fields, but also by (massless) spin-1 fields. [This generalization is formal because, as we 
saw above, spin-1 fields cannot couple exactly to the mass]. If g s denotes the coupling constant of spin-s fields one 
finds @ 

G = EH s .9 s 2 , (80) 



1 E S H S 9 2 S 



7 + 1 = ± ' ys = - < s 2 > (81) 

[We denote here g 2 what was denoted g s in Ref. [p9[]. Note the elegant interpretation of 7 + 1 as being half the 
average squared spin of the mediating fields [the weights being defined by the contributions of the fields to the 2-body 
interaction, including the sign which is negative (repulsion) for s = 1]. When gi = and 90/92 = « 2 , eq. ( |3l| ) yields 
©• 

The quantity (3 — 1 parametrizes the possible deviations from general relativity in the non-linear terms. Its 
expression in a general tensor-multi-scalar theory is 



where 



= 1 a a f3 ab a b 
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(82) 



(3 ab = D a D b In A = D a a b = d a a b - j ba a c 



(83) 



is the second covariant derivative (with respect to the cr-model metric 7 a fc) of the logarithm of the coupling function 



Many observations in the solar system have been used to study the post-Newtonian effects present in eqs. ( |70|) 
(78). [A famous example is the secular advance of the perihelion of Mercury, already discussed in Sec. ffl. At present, 



two sorts of experiments stand out as giving the tightest constraints on 7 and (3. [See [|[ [?]] for a discussion of the 
other tests of post-Newtonian gravity]. Time-delay measurements |57j based on the Viking ranging data to Mars 
|58|l , and Very-Long-Baseline-Interferometry measurements of the deflection of radio waves by the Sun |w| , |6(J , have 
allowed one to measure 7, (nearly) independently from (3. In the former case, this is done by considering the time- 
of-flight of an electromagnetic (radar) signal sent from the Earth, actively reflected on a Viking lander on the surface 
of Mars, and received back on Earth. Writing from eqs. fifty f76|) the curved-space equation for the light cone, 
= ds 2 = g^ydx^dx" , one finds that the coordinate time of flight (ar = ct) is given with sufficient accuracy by 



cAt 



1 + (1 + 7)5 



|dx| , (84) 



where V ~ GMq/t. Using an accurate ephemeris for predicting the coordinate positions of the Earth and Mars at each 
coordinate time t, and the transformation between At and the proper time measured by Earth clocks, one can measure 
the coefficient 1 + 7 of the relativistic (or "Shapiro" [57]) time delay J V\dx\/c 3 appearing in eq. (|8^) by analyzing the 
data where the electromagnetic signals pass near the Sun. [The specific time signature of the "Shapiro" time delay 
allows one to separate it from the many other effects present in the leading "Roemer" time of flight J |dx|/c]. The 
final result of this Viking time-delay experiment is p8[ 

7 = 1.000+0.002. (85) 

The same limit was found in the deflection experiment f39f , while a very recent deflection measurement obtained 
a slightly better limit: 7 = 0.9996 ± 0.0017 fc}. 

The second high-precision test of post-Newtonian gravity comes from an analysis of the laser ranging data to 
the Moon §5l]]. In July 1969, the Apollo 11 mission, besides its spectacular aspect of having landed the first men on 
the Moon, left a panel of corner-cube reflectors on the surface of the Moon. Since August 1969 one has accumulated 
regular measurements of the round trip travel times of laser pulses sent from several stations on the Earth [two at 
present; CERGA, France and McDonald, USA] and bounced off an array of 4 lunar reflectors (the first Apollo 11 
reflector has been completed by two other US-made reflectors - Apollo 14 and 15 - and by two French-made reflectors 
deposited by the Russian Lunakhod 1 and 2 missions. Alas the Lunakhod 1 reflector never sent back any echo). 
Because the Earth and the Moon have non-negligible gravitational binding energies [(E gi ' av /mc 2 )Earth ~ —4.6 x 10 -10 , 
(£ ,grav /mc 2 )Moon — —0.2 x 10~ 10 ], eq. fl69j ) shows that they could fall with a different acceleration towards the Sun. 
Computing d\nG/dtp a from eq. (p4[), and Vip a from eq. ( p39| ) one finds 

E^grav 

( ( J a )solf-gravity = /^g _ _ 3) w _ (86) 

rac 

One recognizes here the effect of the second term on the right-hand-side of eq. (fr2|) . Eq. (|8^) means that the combination 
(4/5—7—3) parametrizes the violation of the universality of free fall happening for self-gravitating bodies in theories that 
differ from general relativity ( "violation of the strong equivalence principle" ) . This effect was discovered by Nordtvedt 
ML who emphasized also that laser ranging to the Moon offered an excellent way of looking for the presence of the term 
( j86"l ) HU . [Note that we are working here under the assumption that there is no violation of the "weak" equivalence 
principle associated with the different compositions of the Earth and the Moon] . Indeed, the differential acceleration 
of the Earth-Moon system in the field of the Sun induces a polarization of the Moon's orbit about the Earth. This 
consequence of a violation of the equivalence principle was, in fact, first pointed out by Newton, see section 6.6 of 
]5^] , and first correctly worked out by Laplace f63|]. For recent theoretical studies of this effect taking into account the 
important mixing with solar tidal distortion see ]64|. |65[| . The most recent analyses of the experimental data yield 

4/? -7- 3 = -0.0005 + 0.0011 (87) 

according to Ref. || and 4/3 - 7 - 3 = -0.0007 ± 0.0010 according to Ref. @. 
Combining ( 35 ) and (871) yields the following value for [3 

(3 = 0.9998 ± 0.0006 . (88) 

To end this section, one should mention the fact that Nordtvedt and Will [Q, ||, || have introduced a more general 
"parametrized post-Newtonian" formalism containing, besides the two parameters 7 and /3, eight other dimcnsionless 
parameters, £, a\, 012, 03, Ci> C2 ; C3> C4> associated with other, a priori conceivable, deviations from general relativity. 



In particular, the two parameters ol\ and a-i are associated with a possible gravitational violation of local Lorentz 
invariance (existence of preferred frames). The original theoretical motivation for considering such preferred-frame 
parameters was the idea that gravity could be mediated in part by a long-range vector field (or by other tensor fields). 
If that were the case one would expect the Universe's global matter distribution to select a preferred rest frame for 
the gravitational interaction. However, we saw above that there were theoretical difficulties in constructing consistent 
field theories that are metric and contain vector or tensor fields. The situation is even worse for the other PPN 
parameters that do not seem to come out of any decent field theory. It is anyway a meaningful phcnomcnological 
question to ask whether all the existing data about solar-system gravity suffice to put significant constraints on all the 
PPN parameters. The answer is yes for most of them ft]. Generally speaking the extra PPN parameters are much 
more tightly constrained than 7 — 1 and j3 — 1 (with the exception of a% which is presently constrained only at the 
level |«i| < 1.7 x 10~ 4 [681, ps|] , Jzo|] and which needs new data to be more tightly constrained |o|, 72 ). 

5.4 Theoretical conclusions about weak-field, metric gravity 

A first conclusion is that general relativity is consistent with all existing tests of weak-field gravity (at the Newtonian 
and post-Newtonian levels). We have argued above that the most natural (and probably the only theoretically 
consistent) metric alternative to Einstein's pure spin-2 theory, is a metric theory where gravity couples exactly to 
mass and is mediated both by one massless spin-2 and one or several, massive or massless, spin-0 fields. Within 
this framework the fraction of the gravitational interaction carried by all the scalar fields is constrained at the 1CP 3 
level or better. Indeed, in the case of finite-range scalars one had the constraint a 2 = a a a a < 1CP 3 from Newtonian 
measurements for ranges between 10 m and 10 km, eq. (^3|) (and tighter constraints for other ranges, except for ranges 
< 1 mm) , while in the case of infinite-range scalars (or with ranges greater than the Earth-Sun distance) the limit ( p5| ) 
yields (when using eq. (|79|)) the same numerical level 

a 2 < 10- 3 (89) 

for the constraint on a possible admixture of spin-0 in the gravitational interaction. 

Note that, while the existing observational limit on the post-Newtonian parameter 7 — 1 yields a rather strong 
constraint on the most natural theoretical alternatives to general relativity, the observational limit (|8^) on the other 
phcnomenologically independent post-Newtonian parameter (3 — 1 gives only a very modest supplementary constraint 
on these alternative theories. Indeed, in the (most favourable) case where a 2 — ^aboPol' is equal to 10~ 3 the limit ( pf ) 
is only telling us that the largest eigenvalue of the matrix (3% — 7 ac ftb, eq. (|83|), must be smaller than about 1.2. 

These conclusions raise several questions. The first one is to know what are the prospects for probing with higher 
precision possible deviations from general relativity in the Newtonian (tests of 1/r 2 law) and post-Newtonian regimes. 
Concerning the Newtonian tests there is certainly room for improvement, coming e.g. from the use of superconducting 
gravity gradiometers that can directly test whether the Laplacian of the gravitational potential vanishes or not |7^] . 

Concerning post-Newtonian tests, the Stanford gyroscope experiment |7^] (now called GPB, for Gravity Probe 
B) aims at measuring the velocity-dependent gravitational effects ( "gravitomagnetism" ) with a precision corresponding 
to the 10 -5 level for 7, and thereby a 2 . A similar precision on 7 will be reached by the successor of HIPPARCOS, the 
cornerstone project of ESA named GAIA. Some dedicated missions (such as the high-precision time-transfer project 
SORT proposed to ESA by C. Veillet et al.) might go down to the 7 — 1 ~ 10~ 7 level. The prospects for improving 
post-Newtonian orbital tests by using artificial satellites are discussed in ref. [Q . An improvement in post-Newtonian 
measurements, below the present level ( |S9"1 ) is important in the light of a recent study of tensor-scalar cosmological 
models [Q. It is found that tensor-scalar metric theories generically contain a natural attractor mechanism tending 
to drive the world toward a minimum of the coupling function A(ip), i.e. toward a state close to a pure general 
relativistic one (a m i n = (9 In A/d(p) m in = 0), with the rcdshift at the beginning of the matter-dominated era providing 
the measure for the present level of deviation from general relativity. The numerical estimates of Ref. indicate 
values of 1 — 7 > 2(0/0. 1)~ 3 / 2 x 10 -5 where fl = j0 mattor / j0 crltlcal j s the ususal dimensionless measure of the average 
mass density in the universe. 

Another natural question is to know whether it will become possible in the future to probe weak-field gravity 
at the next level in the post-Newtonian expansion, i.e. at the (v/c) i level in the Lagrangian ( |70| ) ["Second post- 
Newtonian level"]. Some theoretical studies |76], [7?], [78| have generalized the parametrized post-Newtonian formalism 
to the second post-Newtonian (2PN) level. In particular, the study of tensor-multi-scalar theories shows that there 
appear two and only two new parameters at the 1/c 4 level, namely J78[ 

ft = a a (3 b a (3 c b a c , (90) 
ft = a a a b a c D a (3 bc . (91) 



^From a phenomenological point of view the two free 2PN parameters ft and ft can be considered as new independent 



parameters whose values must be obtained, along with the values of 7 and j3, by fitting the data to a complete second- 
post-Newtonian-accurate model of the solar system. If the best fit gives (7, (3, (32,(3') = (1, 1, 0, 0) that will confirm the 
validity of general relativity at a deeper level than is presently achieved. Ref. j78| has shown that binary pulsar data 
(see below) already give strong constraints on these 2PN parameters: |/?2| < 6 x 10 -3 , \(3'\ < 7 x 10~ 2 . However, from 
a theoretical point of view, the results (^,^l|) tell us that such a deeper test (2PN versus 1PN) is not really probing 
new, independent theoretical possibilities (at least among the most natural alternatives to Einstein's theory). Indeed, 
if we assume that all the scalar fields carry only positive energy, which means that the cr-model metric j a b is positive 
definite, then the limit (j39|), namely "f a bO a a b < 10~ 3 , is severely constraining the magnitude of all the individual 
coupling constants a a . As we expect (3 a b and D a (3b c to be of order unity we are led to conclude already from the 1PN 
experimental results that P2 and (3' must be numerically small compared to one. In other words, the message seems to 
be that it is theoretically more important to put more effort in determining the "old" 1PN parameters, especially 7, 
rather than in trying to look for 2PN deviations from general relativity. [Practically speaking, this means fitting the 
data to a theoretical description of gravity given by the full 2PN limit of general relativity plus the 1PN deviations 
parametrized by 7 — 1 and (3 — 1] . 

6 Testing the strong and radiative gravitational field regimes 

6.1 Binary pulsars as laboratories for probing strong and radiative gravitational fields. 

All tests of relativistic gravity discussed above have been performed within the solar system, i.e. within conditions 
characterized by a slowly changing and extremely weak gravitational field. For instance the relativistic gravitational 
potential of the Sun at the Earth distance GMq/^tq^ « 1CP 8 , which is of the same order of magnitude as the 
velocity-dependent effects cx (t> orbltal / c) 2 ~ [10~ 4 ] 2 . The measured relativistic effects are but small perturbations to 
Newtonian expectations (as is well expressed by the terminology of "post-Newtonian" regime) . In other words, all the 
solar system tests have probed the gravitational interaction only in the combined limit of weak and quasi-stationary 
gravitational fields. Therefore, even when the experimental accuracy is high these tests have an important qualitative 
weakness: they say a priori nothing about how the "correct" theory of gravity might behave when the gravitational 
field is very strong (such as near a neutron star or a black hole) or very rapidly varying (as in radiative phenomena). 

Fortunately the discovery of binary pulsars in 1974 ]79| ] opened up an entirely new testing ground for relativistic 
gravity, giving us for the first time an experimental handle on the regime of strong and/or radiative gravitational 
fields. Pulsars in gravitationally bound binary orbits provide nearly ideal laboratories for the testing of strong-field 
gravity: being neutron stars, they have surface gravitational potentials GM/c 2 R G(1.4M0)/c 2 (lO km) ss 0.2; they 
move with mildly relativistic velocities (v/c ~ 10~ 3 ) through a repetitive cycle well suited to experimental averaging 
techniques; and they emit periodic pulses of radio noise, detectable over interstellar distances, in some cases as stable 
as the ticks of an atomic clock. The many orders of magnitude separating the self-gravitational fields of pulsars 
(GM/c 2 R « 0.2) from that of the Earth (GM m /c 2 R m ~ 10" 9 ) or even the Sun (GM Q /c 2 R Q ~ 10~ 6 ), and their 
closeness to the black hole limit ((GM/c 2 R)bh = 0.5), make it clear that they give us access to strong-field gravity. 
The modest increase in orbital velocity ((v/cjpsR ~ 10~ 3 versus {v/c)q « 10~ 4 ) does not indicate clearly why they 
can also give us a handle on radiative gravitational phenomena. This comes from the fact that the corresponding 
orbital periods -Ppsr are of order of a fraction of a day instead of a year, and that a binary pulsar is made of two 
objects with comparable masses ~ 1.4M Q , while in the solar system the planets are much less massive than the sun. 
Taking these two facts into account one calculates easily that the change in orbital longitude over some given time 
span, due to gravitational radiation damping, is greater in a binary pulsar than in the orbit of the planet Mercury by 
a factor ~ (M /M Mercury ) x (P M crcur y /fpSR) 11/3 ~ 10 16 . 

A last, but not least, advantage of binary pulsar systems over the solar system is their theoretical simplicity 
and purity. On the one hand, the solar system is a very complex dynamical system with many degrees of freedom 
(many more than can be modelled in full detail; in particular the modelling of asteroids is limiting the precision of 
many solar-system relativistic tests). On the other hand, a binary pulsar has essentially only 6 degrees of freedom 
(although the spin degrees of freedom must be accounted for). In many respects, a binary pulsar is the hydrogen atom 
of relativistic gravity, and like its electromagnetic analog it has allowed one to investigate fine and hyper-fine levels of 
structure of the gravitational interaction (including a classical surrogate of the Lamb shift, i.e. radiative effects in the 
orbital motion). 

After the discovery of Hulse and Taylor |7Stl many authors realized the potentialities of binary pulsars for probing 
strong and/or radiative gravitational fields. In the following, we summarize the comprehensive approach of Damour 
and Taylor |8C|] to which we refer for details and references to earlier work. For reviews of the use of pulsars as 
physics laboratories see the special issue of the Philosophical Transactions of the Royal Society celebrating the 25th 
anniversary of the discovery of pulsars |^l| . See also the Nobel lectures of Hulse and Taylor (8^] . 



6.2 Phenomenological analysis of binary pulsar data ("Parametrized Post-Keplerian 
Formalism" ) 

Binary pulsar data consist in the recording of the time of arrival, shape and polarization of successive electromagnetic 
pulses emitted by a pulsar member of a binary system. The binary nature of the system implies that these data contain 
a wealth of information about gravitational physics. Let us first consider the "timing" data (recording of the times of 
arrival of the centers of the pulses). The intrinsic pulse mechanism is believed to be due to the spinning motion of a 
neutron star. Like a rotating beacon atop a lighthouse, the rotation of the magnetosphere structure of a neutron star 
sweeps a radar beam across the sky. If this beam passes over the solar system, one observes from the Earth a radio 
pulse for each turn of the pulsar. The time of arrival on Earth of each pulse must be corrected for the Earth motion 
around the Sun and for the dispersion due to the propagation of the electromagnetic wave in the interstellar plasma. 
After having done these corrections, the sequence of times of arrival gives us a direct handle on the orbital motion 
of the pulsar. If the timing precision is high (say, 1 to 10 microseconds) one can study in detail many aspects of the 
relativistic two-body problem. More precisely it has been shown that all the independent relativistic timing effects 
bigger than or equal to (v° ihltal / c) 2 Pb (where Pb denotes the binary period) can be described by a simple mathematical 
formula common to a wide class of relativistic theories of gravity (the class of boost-invariant metric theories, i.e. the 
tensor-multi-scalar ones). This "timing" formula Q predicts that the time of arrival (corrected for Earth motion and 
dispersion) of the Nth pulse (where N is an integer) reads 



t N -t Q =F[T,{p K };{p PK };{ q PK }} , 
where T is the pulsar proper time (corrected for aberration) corresponding to the iVth turn, i.e. 

N/2n = v p T+- u p T 2 + ± i) p T 3 
(where v p = l/P p is the pulsar frequency), where 

{P K } = {Pb, T , e , ujo, x Q } 

is the set of "Keplerian" parameters, 



(92) 



(93) 



(94) 



{p PK } = {fc, 7 ,P fc ,r, S ,<5 e ,e,i} (95) 

the set of separately measurable "post-Keplerian" parameters [among which 7 denotes a dimensionful time-dilation 
parameter to be distinguished from the post-Newtonian parameter denoted by the same letter], and 

{q PK } = {5 ri A,B,D} (96) 

the set of not separately measurable "post-Keplerian" parameters. The right hand side of eq. (B4) is given by 



F(T) = 
Ar = 
Ae = 
A s = 
A A = 



D-^T + A R {T) + A E (T) + A S (T) + A A (T)\ , 

2; sin a; [cos it — e(l + i5 r )] + x[l — e 2 (l + Sg) 2 } 1 ^ 2 cos a; sinu , 

7 sin it , 

— 2r ln{l — e cosu — s [sin a; (cos it — e) + (1 — e 2 ) 1 / 2 cosw sin it]} , 
j4{sin[a; + j4 c («)] + esinw} + £?{cos[a; + ^4 e (u)] + ecosw} , 



where 



x Q + x(T - T ) , 
e + e(T - T ) , 

and where A e (u) and to are the following functions of u, 
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uj = luq + k A e {u) , 
and u is the function of T defined by solving the Kepler equation 
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(99) 
(100) 
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(102) 
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(104) 
(105) 

(106) 



Although the splitting of F(T) into the various contributions (|98|)-(101) is a coordinate-dependent concept, one can 
loosely say that Ar represents the time of flight across the (relativistic) orbit ("Roemer time delay"), Ae represents 
the combined gravitational and transverse-Doppler redshifts of the pulsar clock ("Einstein delay"), A5 the gravita- 
tional time delay of the electromagnetic signal propagating in the gravitational potential generated by the companion 
("Shapiro time delay"), while A a is associated with aberration effects. 

Each theory of gravity makes specific predictions about how the various parameters (^4|)-(|9^) are related be- 
tween themselves, as well as to the (a priori unknown) masses of the pulsar and its companion. But the essence 
of the phenomenological analysis of pulsar timing data (the so called "parametrized post-Keplerian" formalism) is 
to a priori ignore the existence of such theoretical relations, and to (least-squares) fit the experimental data to the 
formula (j9^) . The net result of this fit will be to extract in a phenomenological manner from binary pulsar data the 
Keplerian parameters (M), together with the 8 independent post-Keplerian parameters (Pq). [For simplicity, we skip 



the discussion of the fate of the parameters (96)] 



This approach has been generalized to the other pulsar data, those concerning the shape and polarization of the 
successive pulses. Namely, one can write analogs of the formula ( p2| ) for the pulsar-phase dependence of the observed 
flux density Sobs(fobs, 0) and linear polarization angle tp(4>). 

SobsKbs, 4>) = G[0 ; {p k } ; {P PK }} , (107) 
m=H[^;{p K };{p PK }} > (108) 

where <f> is the rotational phase of the pulsar, and where 

{p PK } = {A, A, k, k, a, &, ip , k', a', &'} (109) 

is a new set of post-Keplerian parameters, extractable in principle from pulse structure data. 

Summarizing, the parametrized post-Keplerian approach shows that, besides the easily measured Keplerian pa- 



rameters (|94D , up to 19 observable post-Keplerian parameters listed in eqs. \ p3j ) and ( |109| ) can be extracted in a 
phenomenological manner from binary pulsar measurements. Any theory of gravity will predict some specific relations 
linking these post-Keplerian parameters to the Keplerian ones, to the masses mi and of the pulsar and its compan- 
ion, and to the Euler angles A, r\ of the pulsar spin axis. In each theory of gravity we can use 4 of the phenomenological 
observables to deduce the values of mi, m2, A and rj, so that the redundant 15 post-Keplerian observables give us 15 
tests of the relativistic law of gravitation. 



6.3 Theory-space approach to binary pulsar tests: introduction of strong-field param- 
eters ft, /?', /?",. . . ' 
What is the theoretical significance of the 15 possible tests obtained by combining measurements of phenomenological 
parameters ? What are these tests teaching us about gravity, and especially about strong- field and/or radiative aspects 
of gravity ? To answer these questions, it is necessary to generalize to the strong-field regime the alternative-theory 
approach discussed above in the quasi-stationary-weak-field context of the solar system tests. Fortunately, the same 
class of tensor-multi-scalar theories can be used to define a strong-field and radiative contrast to general relativity. 

This was first pointed out in the context of the original Jordan-Fierz-Brans-Dicke theory. There it was shown 
that strong relativistic internal gravitational fields could modify the orbital dynamics already at the "Keplerian" level, 
and generate an a priori strong emission of dipolar scalar waves J84| (with the observable consequence of inducing a 
corresponding orbital period change in a binary pulsar). However, the solar-system tests constrain already so much the 
only free parameter of the Jordan-Fierz-Brans-Dicke theory that this theory no longer provides a sufficient contrast to 
Einstein's theory even for what concerns strong-field induced dipole radiation effects in a system like PSR 1913+16. 
[Recently, it was suggested that the 11-minute binary X-ray source 4U1820-30 could provide a better testing ground 
pq] . However, this system is likely to be perturbed by the gravitational field of the globular cluster in which it resides 
at a level whichprevents one from using the observed orbital period change as a test of gravitation theories; see the 
update in Ref. gj; see also Ref. |^q| ]. 

Recently, the predictions of the most general class of tensor-multi-scalar theories (containing several arbitrary 
functions) have been worked out in detail, with special emphasis on the effects of strong relativistic internal gravitational 
fields on the orbital motion and gravitational radiation reaction in systems of neutron stars ]29] . One of the main 
results of this study has been the finding that there existed (under some assumptions) a strong-field analog of the 
weak-field theory parameters 7 and (3 introduced by Eddington. More precisely, one finds that, when expanding in 
powers of the fractional self-gravity s = — E STav / mc 2 all the strong-field and radiative effects in binary pulsars, the 
coefficients of these expansions depend on an infinite sequence of theory parameters 



7 i,/?i,/3 2 ,/3',/?",/3 3 ,(/3/3'),- 



(110) 



All the parameters ( 110 ) are explicitly calculable in terms of the arbitrary elements entering the action (|4 
instance (the subscript indicating that an expression is evaluated at the present cosmological values tp^) 

7i = (a a 7afca & )o , 

Pi = (a a p ab a b ) Q , 

p2 = {a a (3 b a Pb c a c )o , 

= (a a a b a c D a (3 bc ) , 

p" = (a a a b a c a d D ab p cd ) Q . 



For 

(111) 
(112) 
(113) 
(114) 
(115) 



The first two parameters (111), (112) are equivalent to the weak-field parameters 7 and P [see eqs. ( |79| ) and ([82])]. The 
further parameters P2, P' , P" , p3,- ■ ■ parametrize deeper layers of structure of the relativistic gravitational interaction 
which have been left unprobed by solar system tests. [As was said above, the second layer, /3 2 , P' could be probed by 
solar system tests reaching the second-post-Newtonian level of weak-field gravity; see [}r8[]. In pictorial language, each 
parameter in the list (110) represents an (a priori) independent direction away from general relativity in the space 
of tensor-multi-scalar theories of gravity. The "post-PPN" parameters P2, P' , P", Pz,. ■ .provide a chart for the yet 
essentially unexplored domain of strong-gravitational field effects (both in the motion and the gravitational radiation 
of systems of strongly self-gravitating bodies). 

To give a feeling for the physical significance of the strong-field parameters ( |110| ) let us mention that the La- 
grangian describing the motion of N strongly self-gravitating bodies [at the approximation where one treats exactly the 
strong self-gravitational effects, but works perturbatively in the inter-body gravitational potential GM /c 2 (distance) ~ 
(i> orbltal /c) 2 ] can be written in the same form (|7(i|), ( |74| ) as above if one replaces G by an effective, body-dependent 
gravitational constant Gab, and similarly 7 by jab and P by P BC [A, B,C = 1, . . . , N being body labels; after having 
done these replacements one should discard the term (4/3 — 7 — 3) (E A i£LV /wiac 2 + E g lav /tubc 2 ) which is taken into 
account in the replacement G — > Gab]- The quantities Gab, Jab, Pbc depend on the strength of the self-gravity of 
the bodies A, B and C. For instance, if one expands Gab in powers of the self-gravities of A and B one finds 



Gab 1 1 , v 

—q-= 1 - -r)(c A + c B ) + 



V H ^ h — 7T~ I & 



cacb 



+(-/ + 1)(P - l)(a A + a B ) + 



^1) (32 -Is': - 



(b A + b B ) 



(116) 



In eq. (116), 7 and P denote the usual weak-field parameters, 77 denotes the combination 4/3 — 7 — 3, while ca, a a, 1>a 
denote some "compactness" fac tors of body A which are of order ca = — 2i91nmA/i91nG ~ — 2E A irav /niAC 2 = O(sa), 
cia = 0{s 2 A ), bA = 0{s 2 A ). Eq. (116) shows how the new, strong-field parameters P% and p' appear at order s 2 . Similar 
formulas exist for the self-gravity expansions of jab and P BC - The strong-self-gravity effects in the gravitational 
radiation emission have also been worked out in detail and shown to depend upon Gab , Jab , Pbc an( i a new quantity 
named {cxaPbPc&d)- Numerically, the compactness ca are of order 0.3 for 1.4 Mq neutron stars, to be compared to a 
maximal compactness ca = 1 for black holes. This is a reason for expecting that self-gravity expansions such as (116) 
are useful, even in the strong-self-gravity context of neutron stars. [See however below]. 



6.4 Experimental constraints on strong-field relativistic gravity. 

Among the ~ 700 known pulsars, only the class of (~ 35) "recycled" (millisecond) pulsars furnishes us with potential 
relativistic laboratories. Among the latter, apart from exceptional cases (notably PSR 1855+09 discussed below), only 
the subclass of short-orbital-period, high-eccentricity binary pulsars with neutron star companions provides interesting 
gravitational laboratories. At present the latter subclass contains only two useful systems: PSR 1913+16 discovered 
by Taylor and Hulse in 1974 and PSR 1534+12, discovered by Wolszczan in 1990 @. 

Up to now the phenomcnological analysis of the PSR 1913+16 data has led to the measurement of only 3 
post-Keplerian parameters [Q: k = LoPb/2n linked to the periastron advance [eq. (105)] 7 linked to the gravitational 
redshift of the pulsar clock [eq. (p9|)], and P b [eq. ( |106| )], the secular change of the orbital period. In general relativity, 
these 3 quantities are predicted to be the following functions of the masses mi and m.2 of the pulsar and its companion, 



w OK (mi,m 2 ) 



3n / GMn 



2/3 



(117) 



7 GR (m 1 ,m 2 ) = lx 2 (l + X 2 ) (^.) V3 , (118) 



n 



P GR ( mi ,m 2 ) = - X,X 2 {GMnfl* ^ y/2 . (119) 



M = m 1 + m 2 , Xi = m 1 /M, X 2 = m 2 /M = 1 - X 1} 
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n = 27r/n, P 4 (e) ^ 1+ ^e 2 + ^e 4 - 



where we have denoted 



The theoretical prediction ( |119| ) for the orbital period change comes from studying the secular effets of gravitational 
radiation reaction in a binary system of two strongly self-gravitating bodies [p9| . 

In graphical terms, the simultaneous measurement of the three post-Keplerian parameters di ohs , 7° bs and P° bs 
defines, when interpreted within the framework of general relativity, three curves in the mi, m 2 plane, defined by the 
equations 

w GR (mi,m 2 ) = Co ohs , (120) 
7 GR (m 1 ,m 2 ) = 7 ° bs , (121) 
P b GR (m u m 2 ) = i> bs . (122) 

[When taking into account the finite accuracy of the measurements these curves broaden to three strips in the mass 



plane]. Equations ( 1 1 7)— ( 120 ) thereby yield one test of general relativity, according to whether the three curves meet 
at one point, as they should. As is discussed in detail in [90, ^j], ^2|, general relativity passes this test with complete 
success [at the accuracy level 3.5 x 10~ 3 , given by the (experimental, plus Galaxy-induced ) width of the P b strip]. 

This beautiful success raises at the same time some questions. As P GK is physically due to the radiative structure 
of the general relativistic gravitational interaction, one is certainly entitled to view the Co — 7 — P b test as a convincing 
experimental evidence for the existence of gravitational radiation. However, the rigorous derivations of P GR show that 
the full strong-field structure of general relativity plays also an essential role in determining the simple (weak-field-like) 



formula (119). The same remark applies to the two other formulas (117) and (118). This is precisely because of this 
entangling of various structures of relativistic gravity that it is useful to analyze the Co — 7 — P b PSR 1913+16 test 
within the more general theory-space approach. This analysis shows that, in spite of its impressive accuracy, this test 
can be passed by theories that deviate significantly from general relativity. 

Fortunately, the recently discovered binary pulsar PSR 1534+12 gives us an independent handle on strong-field 
gravity. The phenomenological (parametrized post-Keplerian) analysis of the PSR 1534+12 data allowed one to e xtra ct 



4 independent post-Keplerian parameters: Co, 7, r and s. [The latter two entering the gravitational time delay (100)]. 
Within each theory of gravity Co, j, r and s are predicted to be some specific functions of the two masses mi and m 2 . 
Therefore these 4 phenomenological measurements define 4 curves in the (mi,m 2 ) plane of the masses of PSR 1534+12 
and its companion (beware that this is a different mass plane than the one associated to PSR 1913+16). This means 
that we thereby get 4 — 2 = 2 tests of any theory of gravity, according to whether the four curves meet at one point. 
One finds that general relativity passes these two new tests with complete success. It is important to note that these 
tests concern the quasi-stationary, strong- field regime without mixing of radiative effects. At present, the accuracy of 
these strong-field tests is not very high, but numerical simulations show that they should steadily improve as more data 
become available. The system PSR 1534+12 may offer also the possibility of seeing (for the first time) the relativistic 
spin-precession induced by the gravitational spin-orbit coupling [through a careful monitoring of the secular changes of 
the pulse shape]. Indeed, in this system (contrary to PSR 1913+16) the spin axis is significantly misaligned with the 
orbital angular momentum (by at least 8°) [^0|. Recently, it has also been possible to measure (at the 20% precision 
level) P b in PSR1534+12, with a result in agreement with general relativity [Q. 

Shifting from the phenomenological to the theory-space approach, one can ask to what extent all the existing 
pulsar data constrain the possible relativistic theories of gravity, beyond the solar-system data. This question has been 
recently addressed, using as space of theories a specific two-parameter class of tensor-bi-scalar theories, called T(f3' , (3"). 
This class was introduced in Ref. p9[ to describe the two yet unexplored directions in theory space associated with 
the strong- field parameters (3' and (3" , independently of the already explored directions (i.e. independently from the 
weak-field directions 7 and (3, and from the strong-field dipole radiation effects explored in Refs. |84|, |8f|). 

One made use of (i) 10 years of high-quality timing observations of PSR 1913+16, (ii), one year of similar data for 
PSR 1534+12, and (ii) a previously published constraint on a possible strong-field violation of the strong equivalence 
principle [ESI, based on an interpretation of the data of the "non relativistic" binary pulsar PSR 1855+09 p3|. Each 



set of data selects within the two-dimensional plane of theories (/?', (3") some allowed region [Specifically, the region 
where the \ 2 statistics, eq. (||), is smaller than the level corresponding to a (formal) 90 % confidence level]. In this 
way, one gets the following three regions in the (3', f3" plane (see Fig. 3 in ref. pW). (i) a thin strip, roughly located 
around the parabola /3" = ((3') 2 , corresponding to the single (0.5 % accurate) u> — 7 — Pb 1913+16 test, (ii) a wide 
potato-shaped region corresponding to the two new (low-precision) tu — 7 — r — s 1534+12 tests, and (iii) the vertical 
strip —1.6 < (3 1 < +1.5 corresponding to the e — Pb 1855+09 test. When combining these three independent allowed 
regions in theory space one gets two interesting results: 

(1) the three allowed regions do admit a non empty common intersection, and general relativity [i.e. the point 
(/?' ,13") = (0,0)] lies well inside this intersection region, 

(2) at the 90 % confidence level the theory parameters (3' and j3" are constrained to lie in a thin parabolic 
segment whose projections on the f3'/3" axes are roughly —1.1 < f3' < 1.6, —1 < (3" < 6. 

/From a quantitative point of view, these limits are less impressive than the ones obtained on the weak-field 
parameters 7 and (3 by using solar-system data. However, they represent our first limits on possible strong-field effects 
in the motion and radiation of systems of neutron stars. One should note also that, from a theoretical point of view, the 
T(f3',{3") class of tensor-bi-scalar theories considered in the previous analysis has the unpleasing feature of containing 
ghost (i.e. negative-energy) excitations. [A feature actually shared by all previously considered strong-field alternatives 
to general relativity, except the uninteresting, because already too constrained, Jordan-Fierz-Brans-Dicke theory]. It 
was thought in Ref. j2j| that the presence of ghosts (i.e. the indefiniteness of the er-model metric j a b) was necessary 
to construct a class of theories obeying the tight weak-field limits (^5|) , (^) , and still exhibiting significant strong-field 
departures from general relativity. [This is linked to the discussion above of one's theoretical pessimism concerning 
possible 2PN deviations, given the existing tight 1PN limits]. Actually, it has been recently discovered |)5) that a 
physically fully satisfactory class of tensor-scalar theories (containing only positive-energy excitations, and satisfying 
the weak- field tests) exhibited non-perturbative strong- field effects (showing up when considering the exact, infinite 
series of self-gravity effects) which allowed strong-field departures from general relativity. The extent to which actual 
binary pulsar data, (and cosmological considerations), constrain these appealing strong-field alternatives to general 
relativity is presently being investigated. 



7 Cosmology 
7.1 Introduction 

All the tests considered above have examined the gravitational interaction on scales between 1 mm and a few as- 
tronomical units (1AU~ 1.5 x 10 8 km). See also |9(| for a discussion of astrophysical tests that we do not consider 
here. In principle, the Universe is providing us with plenty of data concerning the behaviour of gravity on large scales. 
However, most of these data cannot be used as clean tests of the law of gravity because of our lack of a priori knowledge 
of the matter distribution, and/or the low accuracy of the data themselves (especially in certain cosmological data). 
For instance, there are well established cases (rotation curves in the outer regions of many spiral galaxies, velocity 
dispersions in some clusters of galaxies) where there is a significant discrepancy between the mass that we can infer 
from the observed light and the mass needed to hold the system in gravitational equilibrium if Newton's law is assumed 
to be valid. This discrepancy needs not indicate that Newton's faw is at fault (see, however, |97|| , [^8|) because there 
may well be a lot of unseen ("dark") matter in these systems. For reviews on this "dark matter" issue see Ref. ]99fl . 
Fortunately, there are a few cases where one can factor out one's ignorance of the real matter distribution and get 
rather direct tests of the validity of general relativity on large scales. This happens in particular in some cases of 
gravitational lensing of distant quasars or galaxies by intermediate galaxies or clusters of galaxies (For a review on 
gravitational lensing see e.g. Ref. [ 100 1 ) . For instance one observes a giant optical arc near the center of the rich 
cluster of galaxies A370. This arc is an optical mirage, coming from the lensing of the light of a distant galaxy by the 
gravitational field of the intermediate cluster. The radius of curvature of the arc, 9 = 26" ± 2", is therefore a quanti- 
tative way of probing the gravitational field of the cluster. On the other hand, the dispersion of the velocities of the 
galaxies making up the cluster, a = 1300 — 1700 km/s, gives us another quantitative probe of the mean gravitational 
field of the cluster. One can verify whether the general relativistic prediction linking 9 to a 2 (and to the redshifts 
of the distant and intermediate galaxies) is satisfied, without having to know in advance the real matter distribution 
in the cluster. One finds that the test is satisfied within a precision of order 30 %. Therefore, within this precision 
one has verified the general rela tivistic action on light and matter of an external gravitational field on a length scale 
~ 100 kiloparsec. See Ref. [101] for a discussion of the use of A370 and other simple cases of gravitational lensing as 
tests of general relativity at large distances. The typical accuracy of these tests is ~ 30%. 

Let us now turn our attention to the constraints on gravitation coming from cosmological data. A first type of 
constraints comes from considering that our universe certainly went through a very hot and dense phase ("hot big 
bang" ) . The observation of the cosmological microwave background [isotropy and black body spectrum] establishes 



that the temperature of the universe once exceeded 4 x 10 3 K (temperature of ionization of hydrogen). Even if one 
doubts the details of primordial nucleosynthesis, it seems pretty certain that most of the helium in the universe was 
formed during the early big bang, when the temperature exceeded 10 9 K. With less confidence, one can use the 
comparison between the big bang computations of the abundances of light elements (essentially He 4 , with traces of H 2 , 
H 3 and L i 7 ) a nd the present observations to set constraints on the gravitation theory ruling the evolution of the early 
universe 1 102 , 103, 104 1 . [One must however keep in mind the large uncertainties brought by the need to extrapolate 
the observed abundances back in time, as well as the possibility of inhomogeneities in the early universe [105, |106| ]. 
Basically, this comparison is constraining the rate of expansion of the universe during nucleosynthesis. If we assume 
that the number of degrees of freedom of the thermalized matter of the hot big bang around 10 10 — 10 9 K is that given 
by the Standard Model (with three light neutrinos), the nucleosynthesis constraint is essentially giving us a constraint 
on the value of the gravitational coupling constant during nucleosynthesis. One obtain limits of order |AG/G| 10% 
(see e.g. 



107 1 and references therein). 



In principle, one can also constrain gravitation theories by combining measurements of the present expansion 
rate Hq and of the age of the universe To. For instance, it has been argued that the product HqTq was greater than 
0.4 at the "95 % C.L." [108 . However, the large systematic uncertainties in the determination of Hq and To make it 



difficult to assess the significance of such a limit. 

At present, general relativity is nicely consistent with all the above cosmological tests. [Note that, from time to 
time, the H T test has been claimed to present a real problem for general relativity]. As an illustration of the use of 
the various cosmological tests to constrain the space of possible theories one may consult Rcfs 



107, 5^ 



An important issue concerns the possible existence of scalar fields having only gravitational-strength couplings 
to matter. For instance, the "gravitational sector" of string theory contains, besides the standard Einsteinian tensor 
field g^v, some gauge- neutral scalar fields with Planck-suppressed couplings. We refer to these fields as moduli (they 
include the model-independent dilaton S ~ e~ 2 *). The moduli fields are massless to all orders of perturbation theory 
and play a central physical role in string theory in that their vacuum expectation values (VEV) determine the coupling 
constants of the theory: notably the string coupling constant ff string = e* (associated to the string loop expansion), 
the (unified) gauge coupling constant g ^ = S + f(T) + ■ ■ ■ [ 109 1 and the gravitational constant (in string units) 
G ~ a' [e 2 * + ■••]. These crucial properties of the moduli suggest that they may play an important role in cosmology 
and in low-energy gravity. 



7.2 Cosmology of moduli fields 

Let us recall the basic equations of homogeneous Friedmann cosmology (in an "Einstein" conformal frame, i.e. with 
standard kinetic terms for g^ v ) 



ds 2 



g^u dx» dx v = -dt 2 + a 2 (t) 
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The source terms p and p satisfy the energy balance law 
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and are given by a sum of zero-pressure (non-relativistic matter), radiative, field-kinetic-energy and potential contri- 
butions. In the case of one modulus with standard kinetic term: 
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When one type of matter contribution dominates, the pressure over density ratio A = pj p is approximately constant 
(A nr = 0, A ra( j = i, Akinetic = +1) Apotontiai = — 1) and the energy density varies with the scale factor as p oc a~ 3( ^ 1+A -' 



This gives, respectively, a 



,-6 



and a° in the non-relativistic-, radiation-, kinetic-energy- and potential-driven 



expansion. The corresponding expansions follow a power-law, a(t) cx 1 3 ( 1 + A > , when A ^ -1. and an exponential law 
when A = — 1 (potential-driven inflation). [In terms of the conformal time i] ~ J dt/a(t), this becomes a(rj) oc rj 1 + 3X .] As 



there is always some thermal radiation around, note that a kinetic-energy driven expansion is unstable and necessarily 
becomes rapidly dominated by radiation. [This was an argument levelled long ago by Grishchuk against the proposal 
of Zeldovich to use the "hard" equation of state p = p to describe the early universe.] 

Let us also recall some of the puzzling features of our present large-scale universe that we would like to explain 
in a natural manner: 

a. The extreme smallness of the vacuum energy (cosmological constant) on any a priori relevant particle-physics 
mass scale: 

Pvac £ ^ Hi ~ 10- 120 m A P ~ 1(T 41 Aq CD (129) 
(where rap = G -1 / 2 denotes the Planck mass). 

b. The fact that our universe has been expanding for a time > 10 10 yr ~ 10 61 nip 1 without either recollapsing or 
becoming dominated by a spatially negative curvature. This implies that early on the space curvature term k/a 2 



in Eq. (124) was negligible compared to the "time curvature" term H 2 = (a/a) 2 . 

c. The extreme homogeneity of the universe over ~ 10 5 causally disconnected regions at the time of last scattering 
of the cosmic microwave background. 

The inflationary scenarios give a physical explanation of the facts b. and c, but always at the price of some 
fine-tuning of parameters. 



Binctruy and Gaillard [11C| were among the first to ask whether moduli could be useful in cosmology. They 
tried to see whether any of the moduli fields could provide a natural candidate for being an "inflaton" , i.e. for driving 
a sufficiently long stage of exponential inflation through the dominance of its potential energy. They did not find any 



natural candidate for the inflaton among the moduli. Later work, notably one by Campbell, Linde and Olive 111|, 
stressed the specific obstacles to a successful inflationary scenario brought by the existence of the dilaton. In particular, 
instead of driving an exponential inflationary expansion, a (string-frame) constant energy density drives the dilaton 
towards large negative values (corresponding to weak couplings), while the universe expands only as a small power 
of time. Another problem linked to the dilaton is the shallowness of the nonperturbative potentials it might acquire 
|112|] . This makes it difficult to see how (without fine-tuning the initial conditions) its potential can fix the VEV of 
the dilaton at a reasonable value. 



On the other hand, some interesting, qualitatively new, features of string cosmology have been explored 112] 



114], 115], ]llq] . In particular, Veneziano and Gasperini [115], [ 1 1 6 1 , motivated by the "scale-factor duality" of the 



tree-level string effective action 

r d 4 x^e- 2 * R(g)+4{V$Y (130) 



(i.e. the symmetry a(t) — + a -1 (t), Q(t) — > — 31na(t)), introduced a "pre-big-bang scenario" in which our present 
stage of decelerated expansion was preceded by a "super-inflationary" stage of accelerated expansion driven by the 
kinetic energy of the dilaton. In the string conformal frame (i.e. using the cr-model metric g^v), this pre-big-bang 
solution reads (t < 0) 

a oc (-t)- 1 ^ , <J> = _i±^l n( _t) , (131) 
while, in the Einstein frame (metric t^ t „ = e~ 2 * g^) it corresponds to an accelerated contraction 

a oc (~t) 1/3 , $ = — 7=ln(-t) . (132) 
v3 

[Note that the fact that it contracts in Einstein units ensures the stability of this kinetic-energy driven dynamics 
against the unavoidable presence of thermal radiation.] This scenario can provide a large amount of inflation and 
several of its possible observational consequences have been discussed: gravitational waves, relic dilatons, generation 



of primordial galactic magnetic fields [117|. The main shortcoming of this scenario is that it postulates, without being 



able to describe, the existence of a strong-curvature transition between the pre-big-bang stage (131) and the standard 
Friedmann-Gamow hot big bang. 

On a less ambitious vein, one can ask whether the presence of moduli, i.e. the existence of scalar fields with 
Planck-scale natural range of variation, can help in solving the endemic fine-tuning problems of potential-driven 
inflationary scenarios. Let us consider the action describing the dynamics of gravity (described in the Einstein frame) 
and an arbitrary number of moduli fields having a potential V(ip a ) [a — 1, . . . , n] 



S 



Jd*x^l^R(g^)-^ 7ab (<p c )g^d^ a d v <p b -V( V a )} . (133) 



Here fhp = mp/y/iir = (inG)^ 1 ^ 2 is a reduced Planck mass, and the fields ip a are dimensionless. Following (43], it is 
useful to combine the Friedmann equations for the scale factor a(t) with the equations of motion of the moduli ip a (t) 
to write an autonomous equation describing the evolution of the ip's in terms of the parameter 
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Hdt 



dt = In a + const 



measuring the number of e-folds of the expansion. This yields (when k = 0) the simple-looking equation 

■< v + 2^ = -V y ln|y( ¥3 )| , 
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(135) 



where ip' = dip /dp and where ip'/ ov denotes the covariant derivative of ip 1 with respect to the cr-model metric 
7ah(<p) dtp a dtp b . 

The generic solution of Eq. ( |135| ) can be easily grasped from a simple mechanical analogy: a particle with 
position ip and velocity-dependent mass m(ip') = 2/(3 — ip 2 ) moves, in p-time, in a curved manifold submitted to the 
external potential ln| V(ip) | and a constant friction —2ip' ' . If the curvature of the effective potential In V((p) is sufficiently 
small, more precisely if, in the one-scalar case, 



d 2 \nV 



6-±{d v \nV) 2 



< 1 



the motion of ip is rapidly friction-dominated: 



2^ ~ -V„ \nV(ip) 
dp 



In the one-scalar case, Eq. (137) directly gives the number of e-folds as a function of the "inflaton" ip 

2 dip 



N = J dp 

Note that the "scale of inflation" mj , such that 



d v \nV 



V(ip) = mjv{cp) 
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(139) 



drops out completely from equations (135 )- ( Jl 38| ) . Only the dimensionless logarithmic shape \nv(ip) matters. This 
shape needs to be sufficiently flat for inflation to continue during N > 65. Here, one finds the first need of a fine 
tuning: d v lnv(ip) <C 1 to ensure 



Banks et al. [118] emphasized, however, that it is relatively favourable to have a modulus field as inflaton, as a 
canonical scalar field <j> varying on a typical range / < rhp (i.e. V(<j>) = mj v((p/f)) would imply a number of e-folds 
smaller by a factor (/ /fh p ) 2 : N = (f /rhp) 2 f 2dip/d ip \nv(ip) [where ip = <j>/f]. 

Even when the inflaton is a modulus, it remains to find a natural explanation for having a large number of e-folds, 
Eq. (140), i.e. for starting the evolution on a very flat region of \nv{ip). Several mechanisms have been proposed to 
this end. For instance: [119 envisaged potentials levelling off to a constant value when the (canonical) inflaton takes 
large values; [ 120 1 (see also [121]) argued that quantum cosmology suggests that universes, spontaneously nucleating 
out of nothing, preferably start at a maximum of V(ip) [because the instanton action | | = 3mp/8V(ip) is minimized 
there]; and [118 invoked the use of stringy domain walls or other topological defects fL22fl as a natural mechanism for 
triggering inflation [123 , [124]. However, even if one starts the evolution at the top of a potential barrier, one needs 
an uncomfortably small curvature there: if v(ip) ~ v max [l 
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It remains then to satisfy the strong constraint that the density fluctuations generated by inflation be smaller 
than about 10 -5 . This constraint reads 
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The other constraint that horizon-wave-length gravitational waves be compatible with the observed is otrop y of the 
cosmic microwave background is hew ~ Hi /rhp ~ (mi /rhp) 2 ^ 10~ 5 , and is generically weaker then Eq. (141)because 



of the necessary flatness of lnv((p). The constraint (141) creates a mass-scale problem: First, it seems to ex clude that 
inflation be directly generated at the string scale (even if m s t r i n g ~ 10 17 GeV < mp — 3.4 x 10 18 GeV 125 ). Second, 
the inflationary mass scale it suggests, mj ~ 10 15 GeV, seems totally disconnected from the preferred SUSY breaking 
scale msuSY ~ \/ m 3/2 m P ~ 10 11 GeV with m 3 / 2 ~ 1 TeV. [Note, however, that mi can be -C 10 15 GeV if cosmic 
strings are invoked to generate the needed initial density fluctuations.] 

Even if one does not try to use t he m o duli as infl at ons, the existence of light gravitationally coupled fields 
resurrects the infamous Polonyi problem 1 126 , [127], [128], [129 . In essence the problem is that the presence of light 
particles with very weak, Planck-suppressed couplings causes cosmological problems, either because they decay late 
and generate too much entropy while failing to reheat the Universe sufficiently to restart nucleosynthesis, or because 
they do not decay and overdominate the Universe through the energy density stored in the oscillations of their zero- 
mode in their potential V{<p) ~ | fn\{4> ~ 4>m) 2 ■ The three crucial parameters at the root of the problem are: (i) the 
initial displacement Acj> of the VEV of <fi away from the minimum <j) m of its low-energy potential; (ii) the mass of 
<j>\ and (iii) the decay rate T of <f>. The problem arises because of the generally expected links: 



AJ>~ m P , T ~ m\/m% 



(142) 



When ( 142 ) hold one finds that the entire range of masses 10 -28 eV ^ ^ 30 TeV is excluded. More precisely, if 
10 -28 eV ;$ rrifp < 100 MeV (so that r < 7J ) the field (f> has not decayed by now and the energy stored in | cj) 2 + | m 2 , (j) 2 
(which decreases oc a~ 3 ) overdominates the Universe, while, if 100 MeV ^ < 30 TeV the field has decayed by now, 
but its decay has reheated the Universe at a temperature Tp < 1 MeV, too small to restart nucleosynthesis, and has 
produced an enormous amount of entropy diluting away the results of any previous nucleosynthesis. The problem 
cannot be evaded by a long period of ordinary inflation as the latter regenerates via long-wave quantum fluctuations 
(which are important for <f> if m$ <C Hi), an unacceptably large VEV for <f> [ 127 1 . The Polonyi problem is a serious 
difficulty for all moduli because, as stressed in [ 130 1 , 1 131 , [ 132 ] , [ 118 1 , current SUSY breaking lore suggests that they (as 
well as their fermionic partners) acquire masses of order m^ ~ OT3/2 ~ 1 TeV, which is uncomfortably below the 30 TeV 
limit mentioned above. [In essence, this mass estimate follows from V{<f>) — Wgugy v(<f)/mp) — my 2 mpv((j)/mp).] 

Some solutions to the modular Polonyi problem have been proposed. In particular, |l33| pointed out that the 
potential V(4>) = mg USY w(</>/mp), which is at the origin of the problematic value m^ ~ m^/ 2 °f m 4>-, might solve 
the problem by generating a brief period of (secondary) inflation (with one of the moduli as inflaton) at a "weak 



scale" expansion rate: Hi ~ V 1 ! 2 jfhp 



SUSY 



/mp ~ ms/2- They find that a few e- folds of inflation with Hi < m^ 



is enough to sufficiently decrease A0, and thereby the energy stored in V{4>). Another type of solution has been 
proposed in [ 1 2 1 1 . This reference (see also 134 ) shows that the cosmological difficulties of the moduli are avoided if 
the mechanism introduced by Damour and Polyakov [Q for fixing the VEVs of the moduli is at work. 

In brief, the point of [Q] was to show that, contrary to what is usually assumed, having one or several of 
the moduli stay exactly massless in the low-energy world can be naturally (i.e without fine-tuning of parameters) 
compatible with existing experimental data. Two conditions must be satisfied for this to happen: 

(i) string-loops effects must generate a non-trivial dependence of the moduli-dependent coupling functions en- 
tering the effective Lagrangian, i.e. the universal multiplicative factor e -2 * entering the tree-level action (3.1) must 
get replaced by various moduli-dependent functions, 
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where Bi(<p a ) = e~ 2 * + /i($,T) admits extrema at finite values of the ip's; 

(ii) there exist some preferred values of the <p's, say ip^, where all the B^s relevant to determining the low-energy 
mass scales (notably the gauge coupling function Bp(ip a ) which determines Aqcd^)) reach (at least approximately) 
an extremum. A simple mechanism for ensuring this property might be the existence of a discrete symmetry (^-duality, 
T-duality) in moduli space. [See [p4| for further suggestions.] 

Under these assumptions, ]44| finds that a "least coupling principle" holds in that the cosmological expansion 
naturally drives the VEVs of the ip's toward ipf}^, where the moduli (classically) decouple from matter. Estimates of 
the small, but non zero, present values of tp — tp m show that they are compatible wit h existing experimental data, 
including the extremely stringent tests of the equivalence principle (~ 10 -12 level). [ 12 1 1 has studied the consequences 
of this mechanism when considering an early stage of inflation. They found it fully compatible with observational 
facts. In particular, quantum fluctuations are inefficient in regenerating a quasi-classical long-wavelength VEV for the 
</?'s. They also pointed out that, if one considers massive moduli, the Polonyi problem is naturally avoided in such a 
scenario because the moduli acquire during inflation an effective mass m v ~ Hi so that they are very efficiently driven 
to the preferred values ip m at which they store no potential energy. Moreover, in this model moduli within a very wide 



range of masses (which, contrary to usual models, include the SUSY-breaking favored ~ 1 TeV value) qualify to define 
a novel type of essentially stable ultra-weakly interacting dark matter. Indeed, one finds that quantum fluctuations 
generate the following contribution to the cosmological closure density in the form of massive moduli: 



If the inflationary mass scale mi ~ (Hi rhp)% ~4x 10 15 GeV, a modulus of mass m v ~ 1 TeV can close the Universe: 
8 Conclusions 

To complete this review one should mention the fact that gravitational wave observations in Earth-bound interferomet- 
ric detectors (LIGO/VIRGO/. . . ) should soon give us access to many new ways of probing the regime of strong and 
rapidly varying gravitational fields. First, simultaneous observations from an array of interferometric detectors can in 
principle verify whether the gravitational waves received on Earth are of the pure massless, helicity-2 type predicted 
by general relativity, or contain other excitations. Second, the amplitude and shape the detected signals will tell us a 
lot about the theory of gravity. For instance, if the correct theory of gravity is a tensor-scalar one presenting the kind 
of nonperturbative strong-field behaviour mentioned above, one expects that the spherically symmetric collapse of the 
core of a star down to a neutron star state will emit strong monopolar scalar waves, with energy flux of order 

J j^spinOwavcs f~l / J \ 2 

G ' d E^A , (145) 



dt c 5 \dt 



where E srav is the gravitational binding energy of the collapsing core. [See eq. |l45| of Ref. [|29|]. The flux (145) is 
expected to be much bigger than the corresponding pure spin 2 result (which depends crucially on the deviations 
from spherical symmetry). However, the detection on Earth of scalar waves will be hampered by a small factor 
a = \/(l — tJ/2 < 0.032 |29|| . Still, an optimist could hope to detect spin gravitational waves from stellar collapses 
(and thereby to falsify Einstein's theory) before being able to detect the spin 2 waves emitted by the inspiralling 
motion of binary neutron stars ! 

The first general conclusion one can draw from the above review of the experimental situation is that Einstein's 
theory of gravity has passed all presently performed tests with complete success. These tests have probed many 
features of the structure of general relativity: the metric nature of the coupling to matter, the slow-motion weak-field 
limit [which, in field-theory language, gives already some limits on the field content of the theory, see eq. (pl[)], the 
effects of strong internal gravitational fields on the orbital dynamics and the gravitational radiation reaction of binary 
systems, and, to a lesser degree, the large-distance and large-time behaviour of the theory. 

Does this mean that one should stop testing Einstein's theory and consider it as definitely proven ? No, if 
one remembers that general relativity has no free parameters (considering that the cosmological constant belongs to 
the "matter" side of Einstein's equations). Any test of Einstein's theory is a potential killer of the theory. A clear 
experimental disproof of Einstein's theory would represent a major crisis for physics. 

It is interesting to discuss whether one can presently think of ways in which a non-general-relativistic theory 
would have naturally passed all existing tests with the same success as Einstein's theory, while still differing from it 
in an essential way. In fact there are several ways in which this could have happened, and we have already quoted 
some. One way relies on the possible existence of short-range contributions to gravity. For instance, the kind of 
(supergravity- motivated) vector partner of g^ v suggested by Scherk |2l| , with gravitational-strength coupling to the 
(effective) masses of the quarks and leptons (gi = V 'AitGnii) would have escaped detection so far if its range A~lm. 
Other ways use the fact that the cosmological evolution of the universe at large could dynamically drive a non-general- 
relativistic theory to a state where its predictions are virtually identical to the general relativistic ones jl3|, Q . The 
study of such models can help us in focussing on certain experiments which are more likely to unravel significant 
deviations from the general relativistic predictions. 

Appendix 

Our signature is — V H — h; we use greek indices to denote spacetime indices (/i, u, ■ ■ ■ = 0, 1, 2, 3) and latin indices 
for spatial ones (i, j, ■ ■ ■ = 1, 2, 3). The flat (Minkowski) metric is denoted / M „ = diag(— 1, +1, +1, +1) = instead 
of the often used rj^. To save writing minus signs we define g = — det(g^). When using general vectorial frames e a 
(not necessarily coordinate ones) we think of the last lower index on the connection coefficients (i.e. S in Tp S ) as being 
the differentiation index: 

V s V a =e 5 (V a )+T? JS V (3 , 



where = e^d/dx v is the 8-th frame vector (V = V a e a ) viewed as a derivative operator. In Cartan language this 
means using the connection one- forms u a p = T a p S 5 where a is the co-frame (< a , e p >= 8f). [See e.g. Ref. |p5| ]. 
Our curvature conventions follow from Cartan's f2 Q ^ = duj a p + aj a a w°g, which yields, when using a coordinate basis 

R% a = d p T» va + T\ p Y\ a -{p^a} 

Then 
and 

R ee g^R^ . 

With these conventions R^ v and i? are positive-definite for the metrics of spheres. 
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